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The quasihole states of several paired states, the Pfaffian, Haldane-Rezayi, 
and 331 states, which under certain conditions may describe electrons at 
filling factor u = 1/2 or 5/2, are studied, analytically and numerically, in 



\ the spherical geometry, for the Hamiltonians for which the ground states are 



known exactly. We also find all the ground states (without quasiparticles) 



■ of these systems in the toroidal geometry. In each case, a complete set of 



linearly-independent functions that are energy eigenstates of zero energy is 
found explicitly. For fixed positions of the quasiholes, the number of linear ly- 
\ independent states is 2"^^ for the Pfaffian, 2^"~^ for the Haldane-Rezayi state; 

these degeneracies are needed if these systems are to possess nonabelian statis- 
■ tics, and they agree with predictions based on conformal field theory. The 

^ ■ dimensions of the spaces of states for each number of quasiholes agree with 

Q . numerical results for moderate system sizes. The effects of tunneling and of 

the Zeeman term are discussed for the 331 and Haldane-Rezayi states, as well 
^ \ as the relation to Laughlin states of electron pairs. A model introduced by 

' Ho, which was supposed to connect the 331 and Pfaffian states, is found to 



have the same degeneracies of zero-energy states as the 331 state, except at 
its Pfaffian point where it is much more highly degenerate than either the 
331 or the Pfaffian. We introduce a modification of the model which has the 
degeneracies of the 331 state everywhere including the Pfaffian point; at the 
latter point, tunneling reduces the degeneracies to those of the Pfaffian state. 
An experimental difference is pointed out between the Laughlin states of elec- 
tron pairs and the other paired states, in the current-voltage response when 
electrons tunnel into the edge. An appendix contains results for the perma- 
nent state, in which the zero modes can be occupied by composite bosons, 
rather than by composite fermions as in the other cases; the system is found 
to have an incipient instability towards a spin-polarized state. 
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I. INTRODUCTION 



Over the past few years there has been renewed interest in fractional quantum Hall effect 
(FQHE) HI states involving pairing at even-denominator filling factors pj-Q. The earliest 
idea |@] was to generalise the Laughlin state by first pairing the electrons into charge-2 
bosons, then forming a Laughlin state of the bosons, for which the filling factor Ub must be of 
the form i/^ = 1/m, m > even. Since the filling factor u of the electrons is related to that of 
the bosons by z/ = , one obtains u either of the form 1/g or 2/g, where in the second case 
q must be odd. For the cases u = 1/q with q even, this produces a fractional quantum Hall 
state at a filling factor not accessible in the usual hierarchy theory |T^. This idea was taken 



up by Haldane and Rezayi 0, using spin-singlet pairs, to produce a candidate to explain 
the observed z/ = 5/2 plateau |TT[] (using the usual notion that filling factors larger than 2 
involve filling the lowest Landau level with electrons of both spin, and then constructing a 
z/ = 1/2 state in the first excited Landau level). The Haldane- Rezayi (HR) wave function has 
a simple structure, and other paired states with analogous structures, for either spin-singlet 
or spin-polarized pairs, were written down in Refs. @J^. In particular, the Pfaffian state 
of Moore and Read is the simplest paired state for spinless or spin-polarized particles. 
The latter authors also argued that the paired states exhibit pairing of composite particles, 
either bosons or fermions, constructed by attaching an odd or even number q, respectively, 
of vortices to the electrons, for filling factor 1/q These objects behave like particles 



in zero magnetic field, and the wave functions of the paired states can be interpreted as 
Bardeen-Cooper-Schrieffer (BCS) paired wave functions, in their position space form. In 
particular, this makes it easy to understand why the HR state is a spin-singlet. It was 
also suggested that there should be low-energy excitations in which composite particles are 
unpaired (but still consist of electrons attached to vortices), as opposed to breaking the 
electron pairs in Halperin's picture. It was further suggested that quasiparticle excitations 
analogous to those of the Laughlin state which in incompressible states correspond to 



vortices in the order parameter , would carry multiples of a half flux quantum, and thus 
charges in multiples of l/2g, rather than 1/g as in the Laughlin states at filling factor 1/q 
(this also results from viewing the excitations as quasiparticles in the Laughlin states of 
charge-2 bosons). Finally, it was proposed that these quasiparticles obey nonabelian 
statistics. In brief, nonabelian statistics requires that there be degenerate states for well- 
separated quasiparticles, and when the quasiparticles are exchanged adiabatically, the effect 
is not merely a Berry phase representing ordinary fractional statistics, but a matrix acting 
within the space of degenerate quasiparticle states. In the present paper, we do not aim to 
exhibit this action directly, but we do aim to show that the quasihole states of the Pfaffian 
and HR states possess the necessary degeneracies, and to give a physical explanation of their 
origin. 

In subsequent work by Greiter, Wen and Wilczek the physics of the formation of 
the Halperin-type paired states, that is Laughlin states of electron pairs, was elaborated, 
using the Moore-Read Pfaffian state as an example, and several of the points made in Ref. 
were repeated. Greiter et al. also introduced a three-body Hamiltonian for which the 
Pfaffian state at z/ = 1 is the exact zero-energy eigenstate. As regards the statistics of 
the quasiparticles, however, Greiter et al. argued, quite reasonably, that Halperin's picture 
would lead to simple abelian statistics of the quasiparticles. While we agree with much of 
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the physical discussion by these authors (including the argument that the Halperin paired 
states will have 4g-fold-degenerate ground states on the torus in the thermodynamic limit), 
we disagree with their use of the three-body Hamiltonian and Pfaffian-based simple wave 
functions to illustrate their points. Other work on this model and even the observation 
by Greiter et al. themselves that there is a six-fold degeneracy of zero-energy states of the 
three-body interaction on the torus, are more consistent with the predictions of Ref. of 
nonabelian statistics and related properties that are connected with conformal field theory 
(CFT) in two spacetime dimensions. For example, there are gapless Majorana fermion 
excitations at an edge of the Pfaffian state 0,^, in addition to the usual charge-fluctuation 
boson excitations, while the Halperin-type state of electron pairs would be expected to 
possess only the latter. The results that will be obtained in the present paper lend further 
support to the belief that the quasiparticle states that are constructed as energy eigenstates 
of the three-body Hamiltonian of Ref. (and its generalizations to be constructed below) 
do possess nonabelian statistics. 

It was also suggested that the nonabelian behavior might be present only at points 
of special symmetry, and not be generic. Clearly, the three-body Hamiltonian might be 
such a point. Although it was argued in Ref. ||^ that nonabelian statistics is a topological 
property that cannot be altered by small perturbations because the ground states involved 
are assumed to have an energy gap for all excitations, this has not been tested. It is clearly 
an important problem, but it lies beyond the scope of the present paper. 

A further development in paired FQHE states was the realization that some members 
(to be referred to here collectively as the 331 state) of another class of states, of which an 
example was introduced by Halperin [Q, also exhibit pairing P JT^JT5|] . These states have 
come under scrutiny because of their relevance for FQHE states in double-layer systems at 
= 1/2 P^JT7|JT^|T^ . They can also be viewed as generahzed hierarchy states []TU|, and 



so are not expected to possess nonabelian statistics, however, they are still distinct from 
the Halperin idea of a Laughlin state of charge-2 bosons. We will discuss these states, and 



especially recent work by Ho further in Sec. 

As we have mentioned above, the main purpose of this paper is to check the expectation, 
based on the CFT ideas of Ref. P, that the quasihole states of the Pfaffian and HR states 
possess degeneracies above and beyond those that would be obtained for ordinary Laughlin 
quasiholes, or their generalization to the Halperin-type paired states, and thus to lay the 
groundwork for a demonstration of nonabelian statistics. We do this by constructing zero- 
energy eigenstates of those Hamiltonians for which the simple form of ground state wave 
function is correct. These systems serve as model examples, each of which we may hope 
is typical of a universality class (in the sense of ref. P), though the study of the effect of 
perturbations lies beyond the scope of this paper. (We should mention that quasielectrons 
are expected to obey nonabelian statistics like those for the quasiholes, but it is always 
much more difficult to obtain energy eigenfunctions, of which the wave functions take a 
nice form, for quasielectrons than for quasiholes, and the energies will not be zero, nor 
degenerate, though presumably the degeneracies would be recovered in the thermodynamic 
limit for well-separated quasielectrons.) Some of the results for the Pfaffian appeared in an 



unpublished earlier work [^] (see also Ref. [Q), but the method employed here in general 
is related to that used for the edge states in Ref. p. The results are in full accord with the 
earlier expectations. It will emerge that the degeneracies of quasihole states of the paired 
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states can be viewed as coming from breaking pairs of composite particles and placing the 
unpaired (composite) particles in certain single-particle states that contribute zero to the 
total energy; these are "zero modes" . 

Throughout this paper we will use the terminology "particles" to refer to the underlying 
charged particles in the lowest Landau level, which could be either fermions (such as elec- 
trons) or bosons, and not to the composite particles. For a given Hamiltonian, we will also 
refer to energy eigenstates that have energy eigenvalue equal to zero simply as zero-energy 
states. 

In this paper we mostly work on the sphere; we will briefly review this formalism, and 



results for the quasiholes of the Laughlin state |]T0|. One uses a uniform radial magnetic 
field with a total of A^<^ flux through the surface, and in the lowest Landau level (LLL) 
each particle has orbital angular momentum The LLL wave functions on a sphere are 

usually written (in a certain gauge |T^) in terms of "spinor" (or "homogeneous") coordinates 
Ui and Vi for each particle i = 1, . . . , N, with Ui = e"^*/^ cos^j/2, Vi = e"**^*/^ sin^j/2 in 
terms of the spherical polar coordinates 9i, on the sphere. Since these imply that Ui, 
Vi are not independent complex numbers, it is often more convenient, and will simplify the 
writing, to use a nonredundant parametrization of the sphere by a single complex variable. 
This is done by stereographic projection, which gives the definition Zi = 2Rvi/ui, where R is 
the radius of the sphere. Single-particle basis states in the LLL then take the form z^/{l + 
\zi\'^/4:R'^Y~^'^'f'^'^, where the angular momentum quantum number is = — m. In 

this form, the rotationally-invariant inner product of single-particle states on the sphere is 
given by multiplying one function by the conjugate of the other and integrating over the 
Zi plane with no other Zj-dependent factors inside the integral. Only single-particle basis 
states with m < correspond to LLL functions on the sphere and can be normalized with 
respect to this inner product (the normalizing factors will not be needed here). [Note that 
when and R —>■ oo with N^/R^ fixed, in which case the sphere becomes effectively fiat, 
the basis functions (for m fixed) tend to 2;™'e~'^'' Z'^, the basis functions in the plane in the 
symmetric gauge 0.] Many-particle states can thus be written as 

^ = ^ J](l + V4/22)-a+^0/2) (1.1) 

i 

and ^ must be a polynomial of degree no higher than A^^ in each Zi. Therefore, in the 
following we need specify only \1/ in order to describe a state. 

The Laughlin ground state and quasihole states are exact, zero-energy states for short- 
range pseudopotential Hamiltonians [ffO] of the general form 



H = T.T. VMP^,{N^-M) (1.2) 

i<j Af=0 

in which Pij{L) is a projection operator onto the subspace in which the total orbital angular 
momentum of the particles i and j is L; in the summation, M can be viewed as the relative 
angular momentum. For the LLL states, close approach of the two or more particles occurs 
only when their total angular momentum is large. The parameters Vm are the pseudopo- 
tentials. Note that for spinless particles, only even M occur for bosons, and only odd M 
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for fermions. We will later generalise the projection operator notation to three-body opera- 
tors Pijk, and the subspace onto which it projects will be specified by the values of further 
quantum numbers of the chosen group of particles, such as total spin S, or specific values 
of the z component of the spin of each particle, etc. Every projection operator is always 
normalized in the conventional way, with = P. 

The Laughlin states are zero-energy states for the pseudopotential Hamiltonian in which 
Vm 7^ for M < g and zero otherwise (in fact, the non-zero Vm are usually taken to be 
positive). The Laughlin- Jastrow wave function is 

^LJ = Y[i^^-Zjy. (1.3) 

i<j 

Clearly q must be even when the particles are bosons, odd when they are fermions. The 
number of fiux is then = q{N — 1), and the filling factor, u = N/N^, tends to 1/q as 
the number of particles —>■ oo. We will always use the integer g > as the parameter 
specifying the filling factor u = 1/q. In this state, any two particles have relative angular 
momentum M > q ||10[, so it is annihilated by H. This property is preserved if the state 
is multiphed by the quasihole factors U{w) = lli{zi — w), which change the fiux by one 
quantum. These factors can be expanded in powers of each w to obtain the elementary 
symmetric polynomials in the Zj's, 
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n<«2<---<«n 



:i.4) 



which are linearly-independent operators, and the states obtained by multiplying in several 
of these factors span the full space of zero-energy states for each number of fiux = 
q{N — 1) + n, where n is the number of quasiholes P,p!0|. This results from the standard 



fact about symmetric polynomials that they can all be obtained as sums of products of the 
elementary symmetric polynomials. The space of states obtained in this way is equivalent 
to that for N bosons on the sphere in the lowest Landau level with n fiux, or n + 1 orbitals. 
This can be viewed as the q = case of Laughlin's states, which applies since the dimension 
of the desired space of states is independent of q. The dimension of the space is therefore a 
binomial coefficient 




Also, the expansion of 



:i.5) 



niK^^-^^) (1-6) 

i k 



in sums of products of symmetric polynomials in the w's shows that, when the w^s are 
viewed as the coordinates of bosons [T^, the space of available states for these bosons. 



which behave as if in their LLL with N as the number of fiux, exactly coincides with the 
space of zero-energy quasihole states. The dimension of this space is then given by the 
formula for n bosons in iV -|- 1 orbitals, which is the same binomial coefficient (|1.5|) . The 
equivalence of these viewpoints is the basic duality between bosonic particles and vortices 
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within the LLL; it is analogous to the particle- hole transformation for fermions. The count 
can also be performed by using the q = 1 (fermion) case instead of the g = (boson) case. 
It then gives the number of states for N fermions in A*"^ + 1 = N + n orbitals, or for n holes 
obeying Fermi statistics in the same number of orbitals, and these are the same number 
( |1.5[ ). We will often refer to the dimension of the spaces of zero-energy states we find in this 
paper simply as the number of zero-energy states. 

We now summarise the contents of the remainder of this paper. In Sees. |I|, |TT|, and 
rv| we study the quasihole states of the Pfaffian, HR, and 331 states on the sphere, for the 
Hamiltonians for which these ground states are exact. We find explicit wave functions for 
all the quasihole states, and count them to exhibit in particular the degeneracy that occurs 
even when the positions of the quasiholes are fixed. For the Pfaffian and HR states, this 
is related to nonabelian statistics, while for the 331 states it results simply from a layer 
quantum number of the quasiholes. The analytical results are confirmed numerically. In 
Sec. [V|, we consider the ground states of the same Hamiltonians on the torus, that is with 
periodic boundary conditions, and obtain the wavefunctions of the zero energy states in all 
cases. In Sec. |V|, we make a modest attempt to discuss the effects of perturbations on the 
states considered, especially the Zeeman term (for systems like HR that include particles of 
both spin) and tunneling (for double layer systems like 331). We make a full analysis of the 



zero-energy states of a model proposed by Ho ^0[, which we show to be compressible and 
thus pathological at the point where the spin-polarized Pfaffian is among the zero-energy 
states. We also rectify this problem by adding further terms to the Hamiltonian. These 
results are again checked numerically. In addition, we mention an experimental test that 
can distinguish the Halperin and other paired states, by using electron tunneling into the 
edge, for example via a point contact. Appendix |^ contains definitions used in Sec. 0, 
and Appendix ^ analyses a further paired state, the permanent state, in which there are 
spin-singlet pairs of spin- 1/2 composite bosons this state is found to be at a transition 
point to a ferromagnetically-ordered state. 



II. QUASIHOLES OF THE PFAFFIAN STATE ON THE SPHERE 

In this section we will obtain all the energy eigenstates at zero energy for the three-body 
Hamiltonian for which the Pfaffian ground state is exact, for arbitrary numbers of added 
flux, that is for any number of quasiholes. The following sections generalize the results to 
the Haldane-Rezayi and 331 states, and (partially) to the torus. 

The Pfaffian state pj, for even particle number A^, is defined by the wave function 

(2.1) 




'3 J ' 



where the Pfaffian is defined by 



1 



N/2 



Pf^^i = 2^/2 (Ar/2)! ^ sgndH M,(2fc^iy(2fc) 

for an N xN antisymmetric matrix whose elements are Mj^; S^- is the group of permutations 
of A^ objects. The filling factor is 1/q. The Pfaffian state is totally antisymmetric for q even, 
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so could describe electrons, while for q odd it describes charged bosons in a high magnetic 
field. For g = 1, it is the zero-energy state of the lowest fiux of the Hamiltonian |^ 

H = V b\z,-z,)b\zi-zu), (2.2) 

%<j<k 

where the sum is over distinct triples of particles. 

For numerical purposes on the sphere, it is more convenient to work with a projection 
operator form of the three-body Hamiltonian, instead of the (5-functions in ( |2.2|) . The closest 
approach of three particles on the sphere corresponds to the state of maximum possible 
total angular momentum for the three. If the particles are bosons, the largest possible total 
angular momentum is (recall that each particle has angular momentum N^/2). Then, 

for the q = 1 case, the Hamiltonian may be taken as proportional to the projection operator 
onto the (unique) multiplet of maximum angular momentum for each triple of bosons: 

^= E VP,,k{3N^/2). (2.3) 

i<j<k 

The same trick works for the three-body interaction of fermions giving the q = 2 case; in this 
case, the maximum total angular momentum of three particles is 3N^/2 — 3. Some numerical 
results for such Hamiltonians were already given in 0. For larger q, these Hamiltonians 
can be generalized, in such a way that the zero-energy states are obtained from those for 
g = 1 by multiplying by Y\{zi — zj^'^ (it is assumed that for q odd, we are discussing 
bosons, and for q even, fermions). The presence of the latter factor implies that they are 
all zero-energy eigenstates of the projection operators for any two particles onto relative 
angular momentum M = 0, 2, . . . , g — 3 (g odd), or M = 1, 3, . . . , g — 3 (g even) [or the 
corresponding total angular momenta N^, — 2, . . . , A*",^ — g + 3, (g odd), etc.]. The space 
of states annihilated by such projections is in one-to-one correspondence with the full space 
of states of the g = 1 case, and the desired three-body projection operator [onto angular 
momentum 3N^/2 — 3{q — 1)] is the unique one that corresponds under this mapping to that 
already mentioned for g = 1. For each g, the Hamiltonian can then be taken to be the sum 
of the three-body and all of these two-body projection operators. A very similar approach 
works for the other Hamiltonians studied in this paper, so that results for higher g can be 
deduced easily from those for the minimal g. These Hamiltonians can also be written in 
terms of delta functions and their derivatives, so as to arrive at a form suitable for use in 
geometries other than the sphere. An attempt at a Hamiltonian appropriate for the Pfaffian 
at g = 2 in the second paper in Ref. ^ is invalid as it annihilates all states. 

In Fig. |T]we show the two-particle correlation function g{r) for the Pfaffian state on the 
sphere with g = 2 for three sizes, A^ = 10, 12, 14. We have plotted the function versus 
the great circle separation r (in units of magnetic length) on the sphere, so that the largest 
possible value of r is half the circumference, and we have normalized the curves in such a 
way that in an infinite system they would approach 1 at infinity. We see that, although 
for A^ = 10 an exponential decay at large distances is not apparent, for A^ = 12 and 14 
the curve appears to be rapidly approching 1 at large separation, and these two curves are 
almost indistinguishable, in the region where both are defined. The correlation length in 
the Pfaffian, which would be defined as the length over which g{r) — 1 decays by a factor of 
e, is apparently quite large. 
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The Pfaffian state is the only zero-energy eigenstate of H at A^^ = q{N — 1) — 1. Zero- 
energy quasihole excitations can be obtained only by increasing A'^^, as for the quasiholes of 
the Laughlin state as discussed in Sec. I, but in this case the basic objects contain a half 
flux quantum each and must be created in pairs. A wave function for two quasiholes was 
proposed in 0; it is 

^Pf + 2qholes(2l, • • • , ^^A^; Wi,W2) = 

1 V- nfil[(^cT(2fc-l) - Wi)(z^(2fc) - 'W2) + (Wi ^ W2)] 

y ^ sgncr- 



X - ^^y- (2-4) 

i<j 

It is clearly the pairing structure built into the ground state which allows insertion of 
Laughlin-like factors 

f{zi, zj; wi, W2) = {zi - wi){zj - W2) + {zi - W2){zj - Wi) (2.5) 

which act only on one member of each pair, and, since the /'s increase the maximum 
angular momentum for each Zi by 1, A^,^ increases by 1. As the quasiholes are, at least 
approximately, located at Wi, W2, they effectively contain a half flux quantum each, unlike 
the usual Laughlin quasihole that corresponds to a full flux quantum. The same structure 
requires that quasiholes are made in pairs, since the wave function must be totally symmetric 
or antisymmetric. When quasiholes coincide, that is, when wi = W2, a. Laughlin quasihole 
is recovered. 

It is clear that by inserting more factors /, with different w's, into the sum over permu- 
tations, a whole host of zero-energy eigenfunctions can be obtained. However, this involves 
dividing the quasihole coordinates wi, . . . , W2n into pairs in an arbitrary way; the resulting 
functions are invariant only under exchanges of the two quasihole coordinates in each of 
these pairs, or under permutatations of the pairs. One must then ask whether all these 
states, of which there are {2n)\/{2"' n\) = {2n — l){2n — 3) ■ ■ ■ = {2n — 1)!!, are linearly 
independent, and also whether all zero-energy eigenfunctions can be obtained in this way. 

For four quasiholes, the three distinct functions obtained from dividing the w's into pairs 
in three distinct ways obey just one linear relation, as we will now show. (These methods and 
results for four quasiholes appeared previously in Ref. [^.) For more than four quasiholes, 
the following method becomes increasingly impractical, and we will instead use a more direct 
method, inspired by the results for edge states of the Pfaffian in §]. 

It is convenient to write the functions in the more general form 

^p{zi, . . . ,Zn]Wi,. . . , W2n) = Pf {^p(,Zi, Zj] 1^1, ... , W2n) / {Zi - Zj)]'^]^^. (2.6) 

Here the $p must be symmetric and of degree n'm zi, Z2 in order to represent 2n quasiholes, 
that is so that = q{N — 1) — 1 -|- n. We could use products of the /'s in eq (p. 51), 
but it is convenient to use the following (these choices are clearly related by taking linear 
combinations). For n = 2, define 
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$1(^1,2:2; Wi, ... ,^4) = (2:1 - Wi){zi - W2){Z2 - W3){Z2 - W4) 

+ {Zi - W3){Zi - W4){Z2 - Wi){z2 - W2), (2.7) 

^2{zi,Z2;wi,. . . ,104) = ^i{zi,Z2;wi,ws, 102,104), (2.8) 
^3{zi,Z2;wi, ...,W4) = ^i{zi,Z2; 101,104,102,103). (2.9) 

The following identity is useful: For any set of complex numbers Oj, i = 1, . . . N, N even, 
> 2, 



Fi{ai-aj) = 0, (2.10) 

since the square of the Pfaffian is a determinant in which any three rows or columns obey a 
linear relation. Set $ = $1 + $2 + '^'sj then using ( p.lO|) it can be shown that 



= $1 - i<l> = ^{Zi - Z2)^[{wi - W4){W2 - W3) + {Wi - W3){W2 - W4)]. (2.11) 

Hence ^2-^3 = - $2), where 

X = {wi- W2){W3 - W4)/{Wi - W4){W2 - W3) (2.12) 

is the cross-ratio. Thus as functions of zi and Z2, $1, $2? ^3 are linearly related. To show 
there are no further relations, consider the limit Wi W2, W3 W4. We find 

$1 ~ {Zi - Wi)^{z2 - 103)"^ + {Zi - W3)'^{Z2 - Wi)^, 
$2 ~ $3 ~ 2{zi - Wi){Zi - W3){Z2 - Wi){z2 - W3), (2.13) 

which are clearly linearly independent. 

For N = 2 particles, it now follows immediately that there are only two linearly in- 
dependent states of the type shown. For arbitrary even N, we still have to prove that 
{^p{zi, Zj;wi,W2,W3,W4)/{zi — Zj)} gives only two linearly independent states for fixed 
w^s. For any p = 1, 2, 3, we use ( |2.11|) and expand the Pfaffian in powers of $p. All terms 
containing more than one factor of $p cancel using the identity (|2.10|) since they contain 



factors {zi — zj). Therefore the N -particle wave functions satisfy the same linear relation as 
the for all A^. A similar argument shows that use of linear combinations of the $p inside 
the Pfaffian leads only to linear combinations of the same states. We note that the linearly- 
independent states can be taken to be the unique state where $p in eq. ( p.6| ) is replaced by 
$, and that where only one factor of $p in the expansion of the wave function is replaced 
by the other $p being replaced by $. The effect of is to cancel the pairing factor 
(-2o-(2fc-i) ~ ^(7(2fc))~^ for the pair on which it acts. Thus there is a "broken pair" in the wave 
function, as in This observation provides the method to generalise these results to any 
number of quasiholes. As the Wi vary, these states span a space of zero-energy 4-quasihole 
states, whose dimension we will find below when we have results for general n. 

We now turn to the method for arbitrary numbers of quasiholes or added fiux quanta. 
We will first write down the functions, then explain why they both span the full vector space 
of zero-energy eigenstates for each n and N, and are linearly independent. The functional 
form was inspired by those in M. The functions are defined as 
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'^mi,...,mp{zi, . . . ,Zn;Wi, . . .,W2n) = 

1 F {N-F)/2 ^, \ 

2iN-FV2^N - F)/2l ^Is; ^Sna n z^,,, [[ _ ^^^^^^^^^ 

xl[{z,-z,y. (2.14) 

i<j 

In this equation, ^{zi, Z2; Wi, . . . , W2n) is defined so as to be symmetric in the tu's, and is a 
generahzation of the function $ used in the n = 2 case: 

1 

$(zi,2;2;ti'l,...,iy2n) = -jr-p^ XI \{i.Zl-'^r{2r-l)){.Z2-Wr{2r))- (2.15) 

Clearly the integers must obey Q < nik < n — 1 for each = 1, . . . , F, since the flux 
is = q{N — 1) — 1 + n, and can be taken to be ordered and distinct, < mi < 1712 < 
■ ■ ■ < mp < n — 1, because of the antisymmetrization by the sum over permutations; thus 
< F < n. Clearly we must also have F < N; for N > n this restriction can be ignored, 
and the analytic formulas below apply in this limit. A similar caveat applies to the other 
paired states below, but will not be explicitly mentioned after this section. These functions 
( p.l4| ) represent pairing but with F fermions left unpaired. One could think of the unpaired 



fermions as resulting from breaking pairs for N even, but the states also make sense for N 
odd (note that — F is always even). These functions are closely analogous to the excited 
quasiparticle states of a BCS paired system, where the unpaired particles usually occupy 
plane waves. A whole spectrum of such excitations is expected also in the paired FQHE 
states for any number of quasiholes (including zero), but these generally have non-zero 
energy. Here we are interested only in the subset that have zero energy for the three-body 
Hamiltonian, which occur only when quasiholes are present. These states contain fermions 
occupying a certain set of single-fermion wave functions 1, z, z"^, . . . , which can be 

viewed as the LLL wave functions for a flux of — 1, or as an angular momentum multiplet 
of angular momentum (n — l)/2. However, the actual spatial distribution of the unpaired 
fermions in these states is hard to calculate since it must take into account the whole many- 
particle wave function. Since the effective magnetic field seen by the fermions is essentially 
zero except in the quasiholes where the density is lower, we expect that the orbitals have 
weight concentrated on the quasiholes. The occupation of these orbitals contributes nothing 
to the energy for our Hamiltonian, so these are "zero modes" . The number of zero modes is 
n, the number of added (real) flux quanta, and not 2n, the number of quasiholes. Thus one 
cannot say that there is a zero mode locally bound to each quasihole. Instead the zero-mode 
wave functions are shared among the quasiholes. 

We should point out that Greiter et al. also stated that in the presence of quasiholes 
there are wave functions with broken pairs. However, the functions they published for 
states with one broken pair, both with and without quasiholes (see eqs. (9) and (10) in 
the first paper in Ref. P and (6.2), (6.8) and (6.9) in the second), vanish identically when 
antisymmetrized. Probably for this reason, the counting of the number of unpaired fermions 
that can be accomodated "naturally", i.e. in zero-energy eigenstates, in the presence of 
quasiholes, is stated incorrectly to be n broken pairs, when the correct answer is n fermions. 
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Note also that the states for unpaired fermions without quasiholes that they give are not 
zero-energy eigenstates (and thus not obviously eigenstates at all). 

Linear independence of the states ( p.l4| ) is easily established, for fixed w's. After remov- 
ing the factor Y\{zi — ZjY, we arbitrarily divide the particles into pairs and let the members 
of each pair approach each other, one after the other, say zi ^ Z2i then ^ ^4, and so on. 
For each limit we examine the leading behavior; clearly the leading behavior for each limit 
may be a single pole, in view of the paired form of the function, or it may be nonsingular. 
If it is a pole, we take the function that multiplies the pole (its residue), which is a function 
of the remaining coordinates to which the procedure has not yet been applied, as well as of 
the coordinates zi = Z2, z^ = z^, . . . , to which it has, and we repeat the process. If the first 
limit is nonsingular, we call the function we started with the zeroth nonsingular residue; if 
the first limit is singular, but the second is not, we call the first residue the first nonsingular 
residue, and so on. Thus the m-th residue, obtained after the m-th limit, may be singular 
or nonsingular in the next limit; if it is nonsingular we can identify the original states as 
having m = {N — F)/2 unbroken pairs. Then the (A^ — -F)/2-th nonsingular residue will, 
by definition, be a nonsingular function of the paired coordinates zi = Z2, z^ = Z4, . . . , 
zn-f-1 = zn-f, the unpaired coordinates Z(^iy^F)+i, ■ ■ ■ , zn, and of wi, . . . , W2n- Since two 
functions, one of which is singular and the other nonsingular in a given limit are linearly 
independent of each other, it follows that states with different numbers F of unpaired parti- 
cles are linearly independent. For states with the same F, we consider the (A^ — F)/2 + 1-th 
residue, as a function of the remaining F coordinates Z(^n-f)+i, • • • , zn- It is just a Slater 
determinant in these variables, and these determinants for distinct sets of are obviously 
linearly independent. This concludes the proof. 

From the wave functions ( 2.14|) , it is straighforward to enumerate the number of states 



that satisfy the conditions, for given positions of the w's. First we note that for quasiholes, 
a state (the original PfafRan ground state) exists for A^ even, but not for A^ odd. For 2n = 2 
quasiholes, there is a unique possibility, both for A^ even (with F = 0) and A^ odd (with 
F = 1 fermion, in the mi = state). For 2n = 4, there are 2 states both for A^ even and 
odd; for A^ even these are the same as the two independent states found before. In general, 
for given n > and F, there are clearly 

(2.16) 

independent states. Summing over the allowed values of F, which are those with the same 
parity as A^, we obtain, whether A^ is even or odd, 2^~^, by a well-known formula for 
binomial coefficients. This number, which is valid for A^ > ra, is exactly the number of 
conformal blocks for 2n spin fields in the Majorana conformal field theory [Q; see also Ref. 




m 



Just as for the quasiholes of the Laughlin state (see Sec. I), there is a (finite) positional 
degeneracy associated with the positions of the quasiholes. The functions for fixed w's are 
analogous to coherent states formed out of the linearly-independent quasihole states. In the 
present case, this degeneracy can be calculated, for a given F and set of m^'s, by expanding 
all the $'s in powers of the w's: 



$(21,22; Wi, . . .,W2n - 



(n!) 



^1^2 + + ^r4"')ei(-w) + . . . + e2n{-w) 



(2.17) 
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Here em{—w) is shorthand for the elementary symmetric polynomials in the w's, with each 
Wi replaced by —wf 

emiw) = Wi^Wi^. . .Wi^. (2.18) 

il<i2<--<i m 

which arise since each w appears at most once in any term resulting from the expansion of $. 
It is known that linear combinations of products of the elementary symmetric polynomials 
Cm, m = 0, . . . , 2n yield all the symmetric polynomials in 2n variables. Thus, when the 
functions in eq. (|2.14|) are expanded in powers of the w's, we obtain all the symmetric 



polynomials in wi, . . . , W2n, in which the degree in any one w is not greater than {N — F)/2. 
The total number of linearly-independent states (as functions of the z^s), for a fixed F and 
fixed set of m^'s, cannot be greater than the number of linearly- independent symmetric 
functions of the w's obtained in this expansion. Notice that, if the w's are regarded as 
the coordinates of some kind of particles, the symmetric polynomials in the w's can be 
interpreted as the states for 2n bosons each of which can occupy one of (A^ — -F)/2 + 1 
orbitals (the orbitals having the single-particle wave functions 1, w, . . . , they 
form an angular momentum multiplet of angular momentum (A^ — F)/A). The number of 
such symmetric polynomials in the w's for the Pfaffian case is thus 



{N -F)/2 + 2n 
2n 



(2.19) 



It is our claim, that for each of these linearly-independent symmetric polynomials we have a 
linearly-independent many-particle state (for each set of m^'s), and so the upper bound just 
obtained is in fact the answer. To establish the truth of this claim, we again make use of the 
residues of the successive limits ^ Z2, 2:3 — > Z4, . . . , as defined above. The nonsingular 
residues have a simple form; they are proportional to 

{N-F)/2 2n F 

n I[i^2^-Wr)Y.^Sr^^I[CiN-F+k)■ (2-20) 

i=l r=l aeSp k=l 

The last factor, involving the unpaired particles, is a Slater determinant; the first factor is 
simply a product of Laughlin-like quasiholes acting on the coordinate of each pair. We have 
thus reduced the analysis to the case of the Laughlin states, where all the quasihole states 
are linearly independent (see Sec. I), and this establishes our claim. 

The total number of linearly-independent quasihole states, fixing only and n, is then 




{N - F)/2 + 2n 
2n 



(2.21) 



(Notice that this expression incorporates both restrictions F < n, F < N.) This number 
is clearly larger than the number that would be expected if the quasiholes behaved like the 
quasiholes of the Laughlin states. In the present case the expectation, based on assuming 
abelian fractional statistics as for the Laughlin states, would be (in view of the half flux in 
each quasihole) that the number would be given by the formula for 2n bosons which may 
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each occupy any of N/2 + 1 orbitals (for even). We can compare our result with this 
number, which is 



/ N/2 + 2n 
\ 2n 



(2.22) 



For n fixed and N tending to infinity, the ratio tends to 2"^^. Again, this represents the 
degeneracy necessary for nonabehan statistics. 

We now give arguments that the states found above are a complete set of zero energy 
states. The general construction of zero-energy states for the three-body Hamiltonian for 
g = 1, or its generalizations for g > 1, was given in the first part of Appendix A of Ref. [§. 
It shows that without loss of generality, zero-energy states are linear combinations of the 
forms (from which we have removed the ubiquitous factor Y{{zi — ZjY) 

^ Sre5jv/2 nS /fc(^^(2r(fc)-l), ^a(2r(fc))) 

sgna- '- ^ -, (2.23) 

where the fk are symmetric polynomials in two variables. For A^ odd we can write a similar 
form with = 1, . . . , (A^— 1)/2, and include for the unpaired particle an arbitrary polynomial 
factor /o(^cr(Ar))- In order to represent states with 2n quasiholes, the number of flux added to 
the Pfaffian ground state must be n, and so the fk (and /o for A^ odd) must be of degree at 
most n in each coordinate Zi ; these symmetric polynomials must then be linear combinations 
of the linearly-independent forms 

^^^42+42^^1^ (2.24) 

in which Q < rii < n2 < n. These clearly span a vector space of dimension |(n + l)(n + 2). As 
in Ref. if any fk{zi, Z2) vanishes at zi = Z2, then it must contain a factor {zi — ,22)^? and 
this pair is broken and will contribute to the unpaired-fermion part of the functions ( p. 141) . 



The subspace of symmetric polynomials for which this occurs is spanned by {zi — -22)^ times 
those in eq. ( p.24| ), here with < rii < n2 < n—2; this subspace has dimension |n(n— 1). The 



quotient of these two spaces, which represents the symmetric polynomials in two variables 
which do not vanish at Zi = Z2, therefore has dimension 2n -|- 1. But we have already found 
a set of such functions while expanding $, eq. (p.l7|) , and there are exactly 2n -f 1 terms in 



this series, which form a linearly-independent set for the required symmetric polynomials. 
Therefore we may now argue that, in the general form (|2.23| ), we may choose the fk to be in 



either of the two sets, that is those that vanish at zi = Z2 and those appearing in eq. ( p. 17] ). 
The unpaired fermion part takes the given form, after use of the antisymmetrization among 
the F particles involved. If fk does not vanish at z\ = Z2^ then it is part of the expansion of 
$ for that pair. Thus we get exactly all the terms that occur on expanding eq. ( |2.14| ). This 



shows that the count of states given is correct, and that they can be viewed as arising from 
the states with quasiholes at fixed positions. 

The numbers for the total number of zero-energy states, which are not resolved into 
angular momentum multiplets, are convenient for comparison with numbers of zero-energy 
eigenstates obtained numerically (using the same three-body Hamiltonian (|2.3|) ), when these 
are summed over all L and L^. We can also work out the decomposition of the functions 
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( p.l4| ) into angular momentum multiplets, using the wave functions ( |2.14|) and applying the 
Clebsch-Gordan decomposition to the angular momenta of the bosons that represent the 
quasiholes in the $ factors, and the unpaired fermions. These numbers, and the angular 
momentum decomposition, are in perfect agreement for moderate sizes. The cases that have 
been checked numerically are shown in Table |. In the Table, we have shown only the results 
for even; similar results, in complete agreement with the analytic formulas, are also found 
for odd. 

As an aside, it is interesting that in this case we have established the linear independence 
of the functions in general, whereas in Ref. we were forced to resort to a case-by-case anal- 
ysis. In fact, the results given here now suffice to complete the proof of linear independence 
of the edge states of the Pfaffian state on a disk. We first note that if we place many (or all) 
quasiholes at the same place by setting all the w's equal, then there are no particles in that 
region of the sphere, and the fiuid has an edge there. If we take the limit of these states 
as — > oo with A^ fixed, then the sphere becomes the infinite plane, and the particles are 
concentrated in a disk at the origin if the quasiholes are at the position on the sphere that 
is mapped to the point at infinity by the stereographic projection. Thus, the problems of 
finding the zero-energy quasihole and edge states are essentially the same. The general edge 
states, that include charge-fiuctuation excitations at the edge, are obtained by letting some 
of the w's deviate from infinity, and expanding in 1/w's gives states that contain symmetric 
polynomials in the paired-particle coordinates, rather than in all the coordinates as in Ref. 
IP). This is merely another basis for the edge states; the number of states at each angular 
momentum level is easily seen to be the same, for sufficiently large A^. In the limit, the 
complete and linearly-independent zero-energy bulk quasihole states, for the case where all 
w's are nearly equal, yield all the edge states of the disk, and this gives a proof that the 
latter are linearly independent. To obtain results for the cylinder, with two edges, we place 
half the quasiholes at each pole of the sphere and take a similar limit. In the limit, the par- 
ticles occupy a narrow band around the equator of the sphere, and the infinite fiux through 
the system makes it equivalent to the cylinder, if we consider states where the particles are 
close to the equator, which again means the if 's must not deviate far from the poles. Note 
that the fact that the particles are spread very thinly along this band in the limit we have 
taken does not affect the construction or counting of edge states, which is independent of 
the aspect ratio of the occupied region of the cylinder. Also, the operation that shifts charge 
from one edge to the other P] is obtained by removing a quasihole from one pole and placing 
it at the other. 

In Fig. |] we show numerical spectra for the q = 2 three-body Hamiltonian, described 
after ( |2.ij| ) (in which the projection is onto L = — 3), with V = 1/6, for A^ = 12 

electrons and A^^ = 2N — 3 fiux, that is, no quasiholes, and u = 1/2. In Fig. ^ we show the 
same but with A^ = 10 and 2 fiux added, so there are 2n = 4 quasiholes. The zero-energy 
states, of which the degeneracies were given in Table |I|, can be seen at E = 0, as can the set of 
angular momentum values obtained in this case. The figure also shows that all higher energy 
states are separated by a significant gap that we expect will survive in the thermodynamic 
limit, as needed for the arguments for nonabelian statistics. To further substantiate the 
claim that this system is incompressible, we include some finite size scaling results. In Fig. 
^, we show the low- lying excited energy levels above the ground state for A^ = 10 and 12, 
versus k = L/R. It can be seen that the spectra lie almost on top of one another. In Fig. 
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we show the size dependence of the lowest excited state, versus for several values of A^, 
including those in Fig. ^ Although it is not quite clear by inspection that this energy gap 
is converging to a constant as — oo, because of the behavior of the points for the sizes 
A^ = 10, 12, we note that in Fig. ^ the lowest-energy points for A^ = 10 seem to lie on two 
sides of a minimum, while for A^ = 12 one lies near a minimum. Hence we believe that these 
results do converge to a finite gap. In Fig. |^, we show the ground state energy for systems 
with = 2{N — 1) — 2, that is two quasielectrons, where there are no zero-energy states 
for the three-body Hamiltonian. These energies do appear to be converging to a finite gap 
as A^ ^ cxD. This energy determines the slope in the total energy density versus density in 
the thermodynamic limit, on the higher- density side of z/ = 1/2. Since the energy density 
is zero on the lower-density side, this represents a discontinuity in the chemical potential at 
u = 1/2 for this Hamiltonian, showing that the system is incompressible. 

III. QUASIHOLES OF THE HALDANE-REZAYI STATE ON THE SPHERE 

The Haldane-Rezayi (HR) state can be written in terms of the coordinates of A^/2 
up-spin particles at 2; j, ... , and A^/2 down-spin particles at 

E sgnff ^ _ I -1 — ^n(^i-^i)'- (3-1) 

Here g > 2 is even for electrons, and odd for bosons, and the filling factor is 1/q. The first 
factor is of course just a determinant. The product over Zi's with no spin labels attached is 
over all particles. The fact that this describes a singlet is discussed carefully in 0. In |^ 
it was pointed out that this state can be regarded as a BCS-type condensate of spin-singlet 
pairs of spin-1/2 composite fermions that consist of a particle and q vortices, from which 
the spin-singlet property can be more easily understood. Some further discussion related 
to the edge states is contained in [||. The HR state is the unique zero-energy state at 
= q{N — 1) — 2 flux of a "hollow-core" pseudopotential Hamiltonian that gives any two 
particles a nonzero-energy when their relative angular momentum is either g — lor<g — 3 
0, again with the assumption that the particles are bosons for q odd, fermions for q even, 
and that q > 2. 

In exact analogy with the Pfaffian state, the wave function for two quasiholes is: 
^hr(4,- 



7V/21 

mi 



Uk= i[{4 - ^i)(4(fc) - W2) + {wi ^ W2)] 

t&Sn/2 (4 ^ 4(l)y ' " ~ ^a(N/2))~ i<j 



Due to the spin-independence of the newly-inserted factors acting on each pair inside the 
sum over permutations, the state is still a spin-singlet and this suggests that the quasiholes 
carry no spin. The two-quasihole state is again a zero-energy eigenstate of the hollow-core 
Hamiltonian. To see this fact, expand the inserted factors for each pair in terms of powers 
of zl ± zlf^k)- -D^^ symmetry between zl and z^^^-, in each factor, it is easy to see that 
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zl — zl(^^-^ must occur to an even power. Thus in the complete wave function, the absence 

of {zl — zlf^i^y^^ for any k, I, and hence the zero-energy property of the ground state, is 
preserved in the quasihole states. 

It is possible to write down directly the forms of all the zero-energy states of the hollow- 
core Hamiltonian, in analogy with those for the Pfaffian and those in Ref. |Q. In terms of the 
coordinates of A^^^ up particles, down particles, the wave functions are linear combinations 
of 

E sgnasgnpn(4(fc)) 'IK^im) ' 11 



TJ- ^g5jy^ fc=l l=\ r=l \Z„{^p^j^r) -^piFi+r)) 

xYli^i-ZjY. (3.3) 

i<j 

Here — = Ni — Fi is the number of unbroken pairs, and we may assume the n^'s, 
mfc's are strictly increasing, as for those in the Pfaffian quasihole states, with < < r;,2 < 
■ ■ ■ < np^ < n — 2, < nil < ^2 < ■ ■ ■ < rrip^ < n — 2; consequently, < F^. < n — 1 for 
each cr =t or |. As written, these states do not have definite spin, but eigenstates of and 
of 5*2 can be constructed. Since the paired particles form singlets, the spin is determined 
by the spin 1/2 unpaired fermions in the sums over a and p, which behave identically to 
ordinary spin 1/2 fermions. Hence the possible spin states are determined by adding the 
spins of particles in different orbitals (labelled by or rrik), with the only constraint that 
an orbital occupied with both an up and a down fermion must form a singlet. The total 
spin 5* of the zero-energy states therefore obeys S < {n — l)/2. Notice that the number of 
fiux in these states is = q{N — 1) — 2 + n. Arguments that the states given are both 
complete and linearly independent, so that our count of states is correct, can be constructed 
straightforwardly as a combination of those in Appendix A of Ref. and in Sec. ^ above; 
we omit the details. 

We may now count the number of linearly-independent states as for the Pfaffian. We see 
that for n = 0, 1, we must have F^ = Fi = 0, and such states exist only for N = N-^ + 
even; they are the ground and two-quasihole states written down above. For n = 2, we find 
two possibilities for both even or odd, like the Pfaffian case. For the general case, we can 
write the number of states for fixed n, F^, F^, and w's in two ways. One is 




(3.4) 

If we sum over F| with F = F-^ + F^ fixed, we obtain a second form, 

(3.5) 




which is the number of states for F fermions in 2(ra — 1) orbitals. The sum over F satisfying 
(— 1)-^ = (—1)^ then yields a total of 2^("~^)~^ = 2^"~^ zero-energy quasihole states of all 
spins for fixed A^ (either even or odd, except for small n as already shown) and w's. 
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We now find the dimension of the space of quasihole states at all spins for fixed A^, 
including the positional degeneracy due to the w's as for the Pfaffian. Then the total 
number of linearly-independent quasihole states, fixing only and n, is 

(3.6) 

Again the ratio of this number to that for positional degeneracy only is 22"-3 as ^ oo. In 
this case, the factor 2^" might give the impression that there is a factor 2 for each quasihole, 
perhaps because each carries spin 1/2. But the result is in fact 2^"~^, which indicates the 
connection with nonabelian statistics. There are only n — 1 zero modes available, similarly 
to the Pfaffian, which can be occupied with either spin, with a final condition on the parity 
of the number of unpaired fermions. 

The two-particle correlation functions for the HR ground state with q = 2 have been 
published previously 0, for 6 particles; they suggest that the correlation length is quite 
large in this sytem also. In Table ^ we show results obtained numerically for zero-energy 
quasihole states of the hollow-core Hamiltonian for q = 2, which agree exactly with the 
general formula, as do the orbital and spin angular momenta. In Fig. ^, we show the 
spectrum of the HR state at z/ = 1/2, when no quasiholes are present, for = 8 particles. 
In the spectra, we have chosen Vi = 1. There does appear to be a gap above the ground 
state, and the form of the low-lying excitation spectrum is similar to that of the Pfaffian. 
In Fig. ^ we show the spectrum when 8 quasiholes are present, for A^ = 6 and A^^ = 12. 
However, we note that because the particles carry spin in this system, it is harder to get 
to sufficiently large sizes to ensure that the results are converging to the thermodynamic 
limit, and larger sizes may be needed in order to prove that the energy gaps are approaching 
constants in this limit. 



IV. QUASIHOLES OF THE 331 STATE ON THE SPHERE 



It is of interest to perform a similar calculation for the quasiholes of the 331 state, even 
though the excitations of this state are known to have abelian statistics. The relation to 
paired states has been discussed in p!^JT3|pil| J^ . 

The 331 state is just one of a family of two-component states, the so-called mm'n states, 
first introduced by Halperin 0, and studied further in Refs. |[IO|Jl^ - |l8| , P^ • Using notation 
I, I for the two components, these states can be written: 



^ / (zl 

^ mm'n V^l 5 



5 ^TVt 5 ^1 5 



T 



i 



n(4-4rn(^i-^/r'n(4-4) 



(4.1) 



k<l 



The general mm'n state (for some values of A^^^, A^^^) is the unique lowest total-angular- 
momentum ground state of a spin-dependent pseudopotential Hamiltonian, that generalizes 
]2|) to the two-component case |[T0[| , 



Kj 



m— 1 m' — 1 

VmP^,{N^ - M, TT) + E KlP^AN^ - M, U) 

M=0 M=0 



n-1 



ij M=0 

(4.2) 



17 



in which the projection operators Pij{L,aa') project onto the spin states cr, a' for particles 
i, j, respectively, as well as onto total orbital angular momentum L. Thus this Hamilto- 
nian gives positive energy to any state in which two T o^' i particles have relative angular 
momentum less than m or m', respectively, or in which an | and a [ particle have relative 
angular momentum less than n. 

For the case when the exponents in these states are of the form m = m' = n = q — 1, 
g > 1, (which give filling factor v = 1/q, and the partial filling factors for f, J, are both l/2g; 
for brevity, we will continue to refer to this class of states with general q as the 331 state), 
then use of the Cauchy determinant identity, 

n(4 - 4) iK^i - n(4 - zir' = det (-^] , (4.3) 

i<j k<l rs Zj j 

allows the ground states to be written in a paired form (i.e., as a spin-indepenedent Laughlin- 



Jastrow factor times a pairing function), similar to the Pfaffian and HR states PJTJ]. In 
terms of BCS-type pairing, this function describes p-type spin-triplet pairing, with each pair 
in the Sz = Q state of a spin triplet |r5| , p0| . 

We will write the quasihole states immediately in terms of broken pairs, (although a 
simple form using Laughlin quasihole operators acting on either t or i spins will be described 
later) : 

\)- cT(iSN-^ k=l 1=1 r=l ^(7(F|+r) >(Fi+r) 

X n(^. - (4.4) 

i<j 

which is particularly similar to the HR case, except that here and mi < n — 1, and so 
< Frj < n. The flux in these states is = q{N — 1) — 1 + n, as for the Pfaffian. 

For the count of states at fixed w's we obtain, again summing over with F = F-^ + 
fixed, 

(4.5) 

which is the number of states for F fermions in 2n orbitals. The sum over F satisfying 
(— 1)-^ = (—1)^ then yields a total of 2^"~^ zero-energy quasihole states of all spins for fixed 
and w's. This is valid for all n except n = 0, in which case there is no zero-energy state 
for odd. 

The result 2^"~^ may also be understood by viewing it as coming from the choice of layer 
index on the Laughlin quasihole operator. Thus quasihole states of zero energy can clearly 
be obtained by multiplying factors 11(^1^ ~ w'^), that act on either particles of spin a =t 
only, or on spin | only, into the 331 ground state. Further, the numbers A^^^, need not 
be equal. However, the fiux seen by both | and | particles must be equal. This leads to the 
condition 2(A''| — N^) = rii — n^, where are the numbers of quasiholes of the two types, 
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in which the factors act on the particles of spin a. For fixed N = + Ni, there are 2^"~^ 
ways to choose the spins of the quasiholes consistent with these conditions, for n > 1. We 
will refer to this construction of the zero-energy states as bosonic, because, as explained in 
Ref. [§, the relation of the two approaches is related to bosonization of fermi systems. 

Returning to the paired, or fermionic, description, the total number of states resulting 
from the positional degeneracy is in this case 




(iV-F)/2 + 2n 
2n 



(4.6) 




This is larger than for the Pfaffian or HR states. It should also be possible to obtain 
this formula from the bosonic description, although we have not done so. The positional 
degeneracy for the quasiholes of each spin is in that approach 

(4.7) 

and the summations are constrained by the fact that for the zero-energy states depends 
on the 

For the 331 state, it is again not too difficult to extend the arguments for the Pfaffian 
or HR cases, to show the completeness and linear independence of the zero-energy states 
found, in the pairing or fermionic form. This is also quite clear in the bosonic form of the 
states. 



V. GROUND STATES ON THE TORUS 

In this section we consider the zero-energy eigenstates of the special Hamiltonians dis- 
cussed above on a system with periodic boundary conditions (a torus), without quasiholes. 
Although we believe that the states we will give span the complete spaces of such states, we 
will not prove this, but will refer to numerical results for confirmation. 



First we briefly review known results for this geometry [^, to flx notation. In the 
Landau gauge A = —Byk, we will take the magnetic length to be 1, and the system to 
be a parallelogram with sides Li, L2, and periods Li, L2e*° = Lit in the complex plane; 
a is the angle between the sides, and r (with Imr > 0) parametrizes the aspect ratio. As 
usual, there are N^p flux through the surface. Many-particle wave functions in the LLL can 
be written 

^iz,,...,ZN) = fiz,,...,ZN)e-^^y^/\ (5.1) 

where / is a holomorphic function; as a consequence of the boundary conditions on \I' it is 
required to satisfy, for all i, 



f{zi, Z2, ■ ■ ■ , Zi + Li, . . .) 



f{zi, Z2, . . . ,Zi, . . .) 



f{Zi, Z2,...,Zi + LiT, . . .) Jcj,^ iTrN42z,/Li+T) 



f{Zi, Z2, . . . , Zi, . . .) 
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Here 0i, 02 represent general twisted boundary conditions, the same for all particles, and 
are the same for all states in the Hilbert space. They can be set to zero without any real 
loss of generality. 

From a general symmetry analysis |2^, any state in the system can be decomposed into 
center of mass and relative motion, as 



.zn) 



Fcra (Z)/rcl(2;i, 



Zn) 



(5.3) 



where /rei is invariant under shifts of all Zi by the same amount, and Fcm is a function of 
Z = Y.iZi only. Given the boundary conditions on /, specified by 0i, 02, there is still some 
freedom of choice in the corresponding phases in the conditions on and /rei, which is 
related to Haldane's k- vector quantum number 
any case, there are always q solutions for Fcm, for filling factor z/ 



and will be useful in the following. In 
p/q {p, q coprime), which 
are mapped into each other by magnetic translations of the center of mass, and consequently 
are degenerate in energy for any translationally-invariant Hamiltonian. 

For the Laughlin state on the torus, the property of vanishing as the gth power as any 
Zi — >■ Zj fixes the relative wave function to be (definitions and results for the theta and 
elliptic functions used in this section are given in Appendix A) 



rei 



i<j 



z,)/L,\Ty^f, 



LJ- 



(5.4) 



The basis states for the center-of-mass wave functions found in Ref. 23 can be rewritten 



(see Appendix A), apart from some constant factors, as p5| 



Fcm(^) = ^ 



£/g + (Ar^-g)/2g + 0i/27rg 
-{N^ - g)/2 - 02/27r 



{qZ/U\qr). 



(5.5) 



Here £ = 0, 1, g — 1 labels the center-of-mass degeneracy, and 0i, 02 have been retained 
for generality. Using the properties given in App. 0, this can be verified to obey the condi- 
tions resulting from ( p. 21 ) that are given in [^. Alternatively one can verify directly that 
eq. ( p.3|) satisfies eq. ( p.2| ). Note that Fcm has q zeroes in the unit cell for Z of sides Li, Lir, 
and linear combinations of these functions span the Hilbert space of a charged particle on 
a torus in a magnetic field with q flux quanta through the torus, in the LLL. The flux seen 
by the particles in these states is A^,^ = qN . 

For the paired states on the torus, the relative motion part /rei must be modified from 
/lj, in particular to change the symmetry under permutations, which can be done at the 
same flux Na^ = qN, by using 



/rei /elliptic /l J • 



(5.6) 



Here /ciiiptic(-2i, • • • , ^a^) is a meromorphic function and, since f^^i must be holomorphic, any 
poles in /eiuptic must be of sufficiently low order that they are rendered nonsingular by the 
zeroes of /lj, which are located on the hyperplanes Zi = Zj. Further, /eiuptic must obey 
2j-iiidependent boundary conditions 



/elliptic (^1; ■ ■ ■ , Zi + Li, . . . , Zn) OC /elliptic (^1 ; 



Zn) 



(5.7) 
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and similarly for Zi ^ Zi + Lit. Thus /ciuptic is analogous to an elliptic function, but in 
complex variables. In general, the phases that are the proportionality factors in (|5.71 ) need 
not be equal to 1, because any phase left by translating /rei can be absorbed (cancelled) by 
modifying the behavior of Fcm, by a shift in 0i or 02 (by vr in the present case) in eq. ( ^.51) 
before setting (pi and 02 to zero. This is the freedom of choice that is related to the k- vector, 
mentioned above. For the paired states, /elliptic is expected to be a periodic generalization 
of the pairing functions, such as the Pfaffian, discussed earlier on the sphere. 

We now review the ground states of the three-body Hamiltonian eq. (|2.2|) and its gener- 
alizations on the torus. These ground states, without quasiholes, were found by Greiter et 
al. M. For the relative motion part we take 



/elliptic, a Pf 





-Zj)/Li\T) 




-Zj)/Li\T) 



(5.8) 



where a = 2, 3, or 4. (The three ratios 'i9a/'i9i(z|r) are essentially the three Jacobi elliptic 
functions sn(2;|r), cn(z|r), and dn(2;|r), up to translations of z and some factors.) These 
then reverse sign under exactly two of the three transformations Zi ^ Zi + Li, Zi ^ Zi + LiT, 
Zi ^ Zi + Li[l -\- r), and are invariant under the third, for any i. The change in sign can be 
absorbed by the effect of a modification of F^ra as explained above, such that the full wave 
function always obeys the same boundary conditions. This structure has the consequence 
that Haldane's k- vector quantum number [^, which lies in a Brillouin zone, and which is 
zero (modulo reciprocal lattice vectors) in all q of the periodic Laughlin states, is nonzero in 
these periodic Pfaffian states. The distinct nonzero values for a = 2, 3, or 4, are determined 
by the behavior under the three translations already mentioned. For the most symmetrical 
choice of r, r = e*'^/^, which gives the system six-fold rotational symmetry in real space and 
in the hexagonal Brillouin zone of k's, these nonzero vectors lie halfway along the shortest 
nonzero reciprocal lattice vectors; only three of these vectors are distinct modulo addition 
of reciprocal lattice vectors. Thus they lie at the midpoints of the edges of the zone, for 
this choice of zone. Thus, all 3q zero-energy states found for u = 1/q have distinct quantum 
numbers and so are linearly-independent (indeed, orthogonal) states. Numerical calculations 
are in agreement with these quantum numbers ||26|| . 

Next we turn to the HR state on the torus. For the HR state at u = 1/2 [q = 2), 
numerical calculations have revealed that there are 10 = 5q ground states of the hollow-core 



model p7| , p8| , and so we expect to find 5q for general q. The full set of wave functions has 
not to our knowledge been obtained previously. 

Following the reasoning for the Pfaffian, we might expect /elliptic to be a determinant 
of 'da'&b/'&fs in zj — Zj, with a, b = 2, 3, 4, or possibly an antisymmetrized combination 
of other functions that each obey the same boundary conditions. (This assumption may 
be too restrictive, since the boundary conditions in fact need only apply to the complete 
function. We will see that it works for HR, but not for 331.) These products of Jacobi 
elliptic functions are not, however, linearly independent. In general, elliptic functions are 
completely determined by the singular part of their behavior near the poles, and by the 
periodicities (see e.g. Ref. ^9]). In the present case, we require that under z ^ z + 1, 
z ^ z + T the elliptic function [with arguments {z\t)] be either periodic or antiperiodic, 

giving four possibilities which we will label ++, H — , — h, , in an obvious notation, and 

we also require that there be a double pole at z = 0, with residue zero. 
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For the ++ case, there is a classic solution to these requirements: the Weierstrass function 
P{z\t). The required functions with the other boundary conditions, which we will denote 
p2, p3, pi, (in the same sequence as before), can be obtained straightforwardly (see App. 

Defining pi = p, four candidates for the relative part of the HR state on the torus are 
obtained: 

/elliptic.a = det pa{{zj - Z^)/Li\t), (5.9) 

and so, on including /lj and Fcm, we obtain 4g states. It is easy to see that they are zero- 
energy states of the hollow-core Hamiltonian. The k values for the cases a = 2, 3, 4 are as 
for the Pfaffian, while a = 1 gives states at k = 0. Since the 4g states have distinct quantum 
numbers, they are linearly independent. 

To obtain the fifth set of q states, we note that, for the ++ case only, the Weierstrass 
function is not the unique solution to the problem posed: we get another solution by adding 
a constant. (For the other cases, this would violate the boundary conditions.) Indeed, we 
could have used a constant in place of pi, but the determinant would then vanish except in 
the case of = 2 particles. If we insert pi + c in place of pi in the determinant, and then 
expand in powers of c, we find that terms of order higher than one in c vanish identically 
because there are rows or columns in the determinant that are equal. The term of first 
order, however, is linearly independent of the zeroth order term (which is /eiiiptic,i) and is 
nonzero. Further, it is a zero-energy eigenstate of the hollow-core model; it clearly has k = 0. 
Explicitly, this function is 

I JV/2-1 

/empties = E sgnasgnp [] pi{{zl^^^^.^ - zI^^^^^)/Li\t). (5.10) 

^^V^ aeS^;, r=l 

We point out that this has an interpretation in terms of unpaired fermions. Unpaired 
fermions must occupy single-particle states that are holomorphic and obey the boundary 
conditions on /eiuptic- For any except ++, there are no such states, and for ++ the only 
state is, again, the constant function. By breaking one pair in /eiuptic,! and putting the two 
fermions (with opposite spin) in this constant state, we obtain /eiuptic.s- We observe that all 
5g states found are spin singlets. The first and fifth sets of q states are linearly independent 
because they have different numbers of poles. The k values found agree with the numerical 
results. 

The existence of five sets of q states for the HR state on the torus is surprising, especially 
as the analysis in Ref. found just 4g sectors of edge states, and one expected 0], on general 
CFT grounds, that the number of sectors (or of primary fields of the chiral algebra of the 
CFT, which comes to the same thing) would be the same as the number of conformal blocks 
in the CFT on the torus. All we can say about this here is that the CFT described in 
Ref. P for the HR state does lead to 4g vacuum sectors on the torus, yet what is actually 
required for constructing a QHE state is a correlator containing many insertions of fields 
in the chiral algebra, corresponding to the particles in the ground state. For the HR case, 
these correlators are the ground states found above, and one set of q sectors has turned out 
to contain two conformal blocks. However one chooses to interpret this in CFT terms (it 
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is probably related to other oddities of the CFT for the HR state |^), the existence of 4q 
sectors in the underlying theory is not in doubt. 

Finally, we turn to the 331 state. Since the pole in /eiuptic is expected to be first order, 
we try 



/elliptic, a — dct 



Mi4 


-4)/Li|r) 




-4)/Li|r) 



(5.11) 



For a = 2, 3, 4, these lead to nonvanishing, zero-energy states for the pseudopotential 
Hamiltonian. We expect, however, that there are 4g states altogether, based on the general 
structure of this abelian quantum Hall state (see, e.g. Ref. |^), and we expect the remaining 
/elliptic to be ++. Clearly the natural choice det di/di is nonvanishing only for N = 2. From 
the theory of elliptic functions in one variable, the constant is the only elliptic function with 
at most one simple pole and these boundary conditions, so we have exhausted the possibilities 
of this structure. However, as pointed out above, /eiuptic need not be an antisymmetrized 
product of elliptic functions in {zi — Zj)/Li that each satisfy the boundary condition; only 
/elliptic itself must have this property, so we should broaden our search. 

The correct solution can be obtained, no doubt, in various ways. One way is to use 
the expectation that /elliptic is a conformal block for a Dirac fermi field on the torus, with 
++ boundary condition in the case of interest (the other functions /eiiiptic.a, a = 2, 3, 4, 
can be viewed in exactly this way). This conformal block is known to exist, and could also 
be obtained by bosonization. We are indebted to Greg Moore who obtained the following 
formula (for = 4) for this function, at our request, by a limiting procedure of first 
considering the conformal block for the fermi field with the general boundary condition, 
that is twisted by e^^'^, e^^^, and taking the limit ipi, "^2 0. The result is (up to factors 
independent of the Zi^s) 

peSpf/2 

where -^[{zIt) = ddi{z\T) / dz. Notice that, like the fifth function for the HR state, there are 
two unpaired fermions of opposite spins, occupying the constant single-particle state that is 
allowed by the boundary conditions. For N = 2, this is all that remains, and this function 
was already noted above. For N > 2 this function contains which is not strictly 

periodic (since no such functions exist for ++ boundary conditions), but obeys 

+ r\T) ~ Uz\t) ^""'^ ^^-^^^ 

and is invariant under z ^ z + 1. When any z] (zj) in /elliptic, 1 is translated by Lit (— Lit), 
the result is /eiiiptic,i(2i5 • • • , P^^^ ^ term that vanishes because the constant — 27ri must 
be antisymmetrized with the constant 1 that represents the missing row and column in the 
determinant. Thus /oiuptic,! is invariant, and we have found the fourth set of q zero-energy 
states. We note that k = for these states, and the quantum numbers of all 4g states are 
distinct, so these states are linearly independent. 
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In this section we have imphcitly assumed that the number of particles N is even. One 
may ask if there are ground states on the torus for N odd also. Such states will have an 
odd number of unpaired fermions. As we have seen, this is possible only in the ++ (or 
k = 0) sector. For the Pfaffian, there is no such sector, so there are no ground states for N 
odd, except for = 1. For 331 the k = states already include two unpaired fermions for 
N even. For odd, there must be just one unpaired fermion, of either up or down spin, 
otherwise the state vanishes. But, on generalizing eq. ( [5. 121 ), one finds that it no longer 
satisfies the boundary condition, except, again, for A^ = 1. Finally, for HR, there is no 
problem constructing a zero-energy state with one unpaired fermion in the k = sector. 
This gives a spin-1/2 doublet of ground states for all odd A^, and we have verified numerically 
that states with these quantum numbers are the only zero-energy states for A^ odd. If these 
ground state wave functions are again interpreted as CFT correlators, then they imply that 
in this -t--f- sector, there are nonzero correlators containing an odd number of insertions of 
the Fermi field. 

Greiter et al. also found formulas for two quasiholes of the Pfaffian state on the torus. 
There should be no great difficulty in extending the results of the present paper to include 
any number of quasiholes on the torus for any of the paired states we have considered. 

VI. EFFECT OF ZEEMAN TERM, TUNNELING, AND OTHER 

PERTURBATIONS 

In this section we address the question of what happens to the degeneracies of the 
quasihole states and of the ground states on the torus when the Hamiltonian is varied from 
the exactly-soluble (for the zero-energy states) forms we have considered up to now. First we 
consider the effect of the Zeeman term —h on the HR state; here /i > is the magnetic 

field gfisB, and az,i is the z component of the spin operator for the ith particle, which can 
be represented by the usual Pauli matrices. The hollow-core Hamiltonian is spin-rotation 
invariant, and its eigenstates are also eigenstates for the total spin. Hence the Zeeman term 
simply has the effect of splitting the multiplets of spin states. For the quasihole states, 
this means that the degeneracy is partially resolved, independent of the locations of the 
quasiholes. The lowest energy states are then those with the largest number of | unpaired 
fermions, and the lowest number of J,. Since F^r < n — 1, and F and A^ have the same parity, 
these states contain F-^ = n ~ 1 and F^ = when A^ and n—1 are of the same parity, and 
either = ra — 1, = 1 or = n — 2, = when A^ and n—1 are of opposite parity. For 
fixed w's these states clearly have a residual degeneracy 1 in the first case, 2(n — 1) in the 
second. (The total spin of these lowest-energy states, which is alternately Sz = {n — l)/2 
or Sz = {n — 2)/2, agrees with the results in Ref. 0].) In the first case in particular, the 
lifting of the degeneracy implies that adiabatic exchange of the quasiholes can produce only 
a phase factor, and so the statistics is abelian. 

We also note here that for the edge states the gapless spectrum of spin-carrying 
fermion excitations at the edge implies that the spin susceptibility of the edge is nonzero, 
and for finite h there will be some t fermions present in the ground state, so there is 
a magnetization at the edge, or one could say the edge is reconstructed in this way; a 
confining potential is necessary to stabilise this effect. If one thinks of the t fermions and 
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their antiparticles, the | fermions, as particles and holes of a chiral Fermi-Dirac sea, then 
the reconstructed state just corresponds to shifting the Fermi energy of the sea. There is 
thus still a gapless branch of fermion excitations for both spin | and | at the edge. Further, 
the degeneracy of 5g on the torus is not split by the Zeeman term, because all the states are 
singlets. 

It is clear that similar effects may be expected for any abelian or nonabelian statistics 
state, when there is a symmetry present, particularly a continuous symmetry: A symmetry- 
breaking perturbation may break the degeneracies and leave some kind of abelian statistics 
behavior. This does not, however, necessarily mean that nonabelian statistics is unstable 
against any perturbation, nor does it mean that the abelian statistics obtained is that of 
some simpler abelian state, such as a Laughlin state of charge-2 bosonic pairs. We also note 
that, for either the Pfaffian or HR state, when there are only two quasiholes, there are no 
degeneracies for fixed quasihole positions on the sphere (so there is nothing to split), yet the 
expectation is that the Berry phase obtained on adiabatically exchanging the two is not that 
in the Laughlin state of charge-2 bosonic pairs ||^, at least if the quasiholes are given by the 
wave functions studied in this paper. In fact, while the breaking of degeneracies does strictly 
mean that nonabelian statistics does not occur in adiabatic transport of quasiparticles, 
other associated properties, such as the spectrum of edge states, are, as we have seen, not 
affected (in general, a splitting of the velocities might occur due to a symmetry-breaking 
perturbation.) Quite similar phenomena are known in the hierarchy states, which have 



abelian statistics, where, for example, SU(A^) symmetry appears in the certain states [jT9 



but the symmetry (or related degeneracies) is presumably broken by the Hamiltonian, both 
for the bulk quasiparticle states and for the edge spectrum. For this reason, we propose that 
such effects do not really represent a change of universality class in the nonabelian systems 
either. 

We now discuss the effect of tunneling between the layers on the 331 states; mathemat- 
ically this is the same as a Zeeman-like term —tJ2i'^x,i (see earlier discussions in [^,^). 
This term can be diagonalized by using a basis of ax eigenstates for each particle, which 
we will label e and o ( for "even", "odd"), given by e = (t + i)/V^, o = (t — D/V^, 
which have eigenvalues +1, —1 respectively under ax,i- In the literature, these states have 
often been denoted "symmetric" and "antisymmetric" , respectively. Unlike the HR case, the 
pseudopotential Hamiltonian for which 331 is exact is not spin-rotation invariant, and the 
energy eigenstates are not eigenstates of = ^J2i<^x,i- Thus the tunneling is a symmetry- 
breaking perturbation, which breaks the conservation of S^- The effect of the tunneling term 
is to modify the states, not merely to split their energies. Nonetheless, when t > is small, 
we may try to use degenerate perturbation theory to understand its effect, which means 
diagonalising the tunneling term in the subspace of the states that have zero energy when 
t = 0; this would give the exact results at first order in t. We are not able to carry this 
out analytically in general, because we do not have the matrix elements of the perturbation 
among these states. We may expect, however, that the degeneracies would be at least par- 
tially lifted, in a similar way as for the HR state with Zeeman, by the following argument. 
In spin space, the pairs in the 331 state take the form + iiti= ^i^j ~ OiOj, using an 
obvious notation for the spinors for the ith, jth particles in a term in which these form a 
pair. Halperin |]15| has argued that the effect of positive tunneling t is to cause a change 
in the 331 ground state, which within a trial-wave-function description causes the amount 
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of 00 in the pairs to decrease. Ho has proposed that this occurs at first order in t. If 
we write the unpaired fermions in our zero-energy states also in the e-o basis, and neglect 
the effect of t on the paired part of the state, its effect would be to split the energies of the 
e and o unpaired fermions, exactly as for the t and | fermions of the HR state, and again 
with the effect of removing most of the degeneracy. For the ground states on the torus, we 
first note that the states with k 7^ are even under layer exchange, which has the effect 
of multiplying by Y\i(Jx,i, while the states at k = are odd (incidentally, this agrees with 
numerical findings [0). The broken pair in those states contains one e and one o fermion. 
Application of the same naive argument as for the quasihole states then suggests that the 
Aq degeneracy is split to 3g by t > 0. However, an accurate calculation should include the 
modification of the rest of the state, and the splitting might disappear in the thermodynamic 
limit. 

If we consider arbitrary small changes in the Hamiltonian from the special forms consid- 
ered in previous sections, then physical arguments like those for the Zeeman and tunneling 
terms suggest that, as the degeneracy arises from breaking pairs and putting the fermions 
into the zero-mode functions, there is no obvious reason why it should not be broken in gen- 
eral. For the bound states in the gap in the vortex cores in conventional BCS superfluids, 
such excited many-particle states have positive energies. However, such arguments may just 
be too naive, because of the modification of the ground state wave function. 

In an interesting paper. Ho ||20| has given a more general interpolation between the 331 
and Pfaffian ground states by varying the spin states of the pairs. Within a paired trial 
wave function description, the effect of tunneling is presumably to reduce the amount of 00 
in the spin part of the wave function of a pair |]15|. If this is done without changing the 
spatial factor {zi — Zj)~^, then when the pairs are purely ee, the ground state is precisely the 
Pfaffian state. Ho claims that this somehow contradicts the "topological" arguments that 
assert that one cannot go continuously between these distinct ground states. However, this is 
really a misstatement; one can always interpolate between any two state vectors in the same 
Hilbert space. The real question is whether the properties of the states, like those considered 
in this paper, can be continuously connected. If one wishes to exhibit a breakdown of the 
"topological" arguments, then it is necessary to show that the interpolation occurs without 
any phase transition, that is without any energy gap for local excitations going to zero 
(which would be a second order transition) or any level crossing of ground states (which 
would be a first order transition). Ho proposed a family of Hamiltonians for each of which 
his corresponding wave function is a zero-energy ground state, but did not show that the 
energy gap is maintained throughout the interpolation. We will now examine this. We 
will show that at the point where Ho's ground state is the Pfaffian, there is an enormous 
degeneracy of other zero-energy states for his Hamiltonian, implying that the energy gap 
collapses to zero at this point, and the system is not incompressible. We will show that, up 
to that point, the degeneracies of the zero-energy eigenstates of his model are the same as 
those of the pseudopotential Hamiltonian for the 331 state, considered earlier. We will then 
consider modifications of his model that remove the pathology. 

Following Ho, we now fix g = 2, so the particles are fermions, and the filling factor is 
z/ = 1/2 (as usual, other cases are similar). To understand Ho's description, first recall that 
for spin 1/2 fermions in the LLL, Fermi statistics implies that if the particles have even 
relative angular momentum, m = 0, 2, . . . , then they must be in an antisymmetric spin 
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state, which can only be a singlet, while if they are in an odd relative angular momentum 
state, m = 1, 3, . . . , then they can only be in a symmetric spin state, which must be a 
triplet. We emphasize that these statements remain true in the presence of any number 
of other particles. If the two particles have total spin 1, then there is a three-dimensional 
complex vector space of spin states, and the spin state can be described exactly by a nonzero 
complex 3-component vector d, of which the magnitude and phase are irrelevant to the 
state. We will not need the detailed definition of d, which can be found in Ho or 



in references on superfluid ^He (see e.g. Ref. 0]). We note only that, for the spin state 
Tiij + ijTj= ejCj — OiOj (as in the pairs in the 331 state), d oc x^, while for CjCj (as in the 
pairs in the Pfaffian), d oc (x^ — ■iXy)/-\/2 (here x^ denotes a unit vector in spin space in the 
z direction, etc.). Notice that the transformation from one state to the other is not simply 
a rotation. Ho's Hamiltonian is a pseudopotential Hamiltonian 

i?Ho(d) = E [VoP^,{N^, S = 0) + ViP,,{N^ - 1, ± d)] (6.1) 

i<j 

which gives positive energy to any two particles with (i) relative angular momentum zero 
and spin zero, or (ii) relative angular momentum one and triplet spin state orthogonal to 
a chosen state specified by d. For the choice of d that corresponds to 331, this reduces to 
the pseudopotential Hamiltonian ([4.2|) used before; pairs may have zero energy and relative 
angular momentum one only if they are in the tijj + J,j|j spin state. For the d vector that 
corresponds to the ee pairs in the Pfaffian, it allows two particles to have relative angular 
momentum one at zero energy if they both have spin e, not if they are eo+oe or oo. However, 
at this d vector, which we call the Pfaffian point, all states of the electrons in which all spins 
are e are zero-energy eigenstates of this Hamiltonian. This follows because there are clearly 
no singlets in such states, so no relative angular momentum zero pairs either. This should 
be no surprise, since no convenient two-body Hamiltonian (like Ho's) giving the Pfaffian 
as ground state and a sensible spectrum is known. Since there are very many spin-aligned 
states at = 2N — 3, Ho's Hamiltonian has a very large ground state degeneracy. There 
might, of course, also be degenerate zero-energy states in which the spins are not all e. We 
will fully analyse this degeneracy below. 

For d gi (x^ — ixy) / ^/2, the degeneracies of Ho's ground state, and of quasihole states on 
the sphere and ground states on the torus, coincide with what was found earlier for the case 
of the 331 state. To construct these zero-energy states for general d, we use a more precise 
notation for the wave functions that includes the spin states. We label the particles i = 1, 
. . . , N, and use |j, |j for spinors which are eigenstates of az^f, the wave function is now in 
a tensor-product space of spatial wave functions and spinors, and must be antisymmetric 
under simultaneous exchange of coordinates and spinors of two particles. For example, the 
ground state on the sphere is 

N/2 /IN 

= E sgn^ n ^^^''"'^"^''^^^ m^^ - ^^r (6-2) 

(which is actually a Pfaffian), where Xi,i(d) is the spin state of two particles i, j that 
corresponds to d, and the product of these factors is actually a tensor product, so that \&d is 
multi-spinor function of the coordinates. This clearly has zero energy. We will call it the Ho 
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state, as it appears in Ref. (in a different notation). In the presence of quasihole factors 
$ in the pairing factors, one may construct zero-energy states with unpaired fermions, as 
for the 331 state, but with the pairs taking the same form as in eq. ( |6.2| ). The unpaired 
fermions can be in either spin state. Even though the ground state (|6.2|) does not have 
definite total Sz, except in the 331 special case, the counting of zero-energy quasihole states 
proceeds just as for the 331 state, and the results are identical to those in Sec. Similar 
results are found for the edge states, which are in one-to-one correspondence with those in 
0, and for ground states on the torus, which are like those in Sec. 0. We expect that these 
zero-energy states are the complete set, except in the Pfaffian limit. This is consistent with 
the continuity of the spectrum as a function of d m a finite size system, which requires that 
the energy levels found in the 331 state must change continuously with the Hamiltonian. 
We expect that the larger degeneracy appears only at the Pfaffian point, because there the 
pairs have the special property of being composed of one spin (e) only, unlike the general 
states This implies that the energy of some excited states decreases as this point 

is approached, so the gap goes to zero. Therefore Ho's Hamiltonian is pathological at the 
Pfaffian point. 

In spite of this pathology at the Pfaffian point, it is still of interest, given that the 
degeneracies of the quasiholes and torus ground states are the same in Ho's model (except 
for the Pfaffian point) as in the 331 state, to ask whether the statistics properties are 
the same. In the 331 state, the structure of these properties is described in terms of a 
U(l) [or S0(2)] quantum number, which in the case of 331 is Sz- Also, excitations with 
opposite values of this quantum number are degenerate, because of the symmetry operation 
of interchanging the layers, which generates a Z2 group. These two symmetries combine 
to make up the symmetry group 0(2), which is the semi-direct product of Z2 and SO (2). 
Ho's Hamiltonian breaks the conservation of Sz, but not the Z2 symmetry, if we consider 
only the family d = x^cos^^ — zxj^sin^, 0<^^<7r/4as proposed by Ho; these correspond 
to the spin state cos(6' — 7r/4)eiej + sin(^^ — 7r/4)ojOj for the pairs, suggested by Halperin 
||15|| , and ^ = is the 331 case, 6' = 7r/4 is the Pfaffian. Consequently, one might think 
the U(l) quantum number is lost. However, the degeneracies of the quasiholes and torus 
ground states are consistent with the presence of this quantum number, which would be 
a "hidden" U(l) symmetry. One might expect this to be, in some sense, the symmetry of 
rotation in spin space about the axis of d, however, since d is complex we must be careful. 
Under rotations of spin space, d is rotated by the action of a (real) orthogonal 3x3 matrix 
in 0(3). In general, there is no rotation that leaves d invariant (up to multiplication by 
a phase), except when d is of certain special forms of which the d vectors for 331 and 
the Pfaffian happen to be examples. Remarkably, even though 0(2) symmetry is broken 
by Ho's Hamiltonian, it seems to be reappearing in the low-energy properties. Indeed, in 
terms of a conformal field theory (CFT) description of the edge states on the cylinder, which 
gives detailed information about the structure of the states [|], the low-energy edge states 
obtained as zero energy eigenstates of Ho's model have the same structure as in for all d. 
If we assume that a CFT description must involve a U(l) theory for charge, together with 
some unitary c = 1 theory, combined by the Z2 orbifold construction of Ref. [§, then the 
theory described in that reference seems to be the only possibility. However, we should be 
cautious about concluding on the basis of these observations that the universality class of 
Ho's model is the same as that of the 331 state, except at the Pfaffian point. We will see 
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below an example in which the full properties of this class do not emerge, and the degeneracy 
can be broken by a perturbation even in the thermodynamic limit. This example emerges 
from further analysis of the Pfaffian point of Ho's model, to which we now turn. 

At the Pfaffian point of Ho's model, in the q = 2 case, it is a two-body projection-operator 
Hamiltonian which gives positive energy to any two particles which either have opposite spin 
and relative angular momentum zero, or if one or both of them is o and they have relative 
angular momentum one. It is known in general how to find the zero-energy states of such 
pseudopotential Hamiltonians; this was already discussed at the beginning of Sec. |^ To 
be zero-energy eigenstates, wave functions must contain the mm'n wave function, eq. ( [4.1| ), 
as a factor, in which for the particular class (but for general q) considered here, we have 
m = g — 1, m' = g + 1, n = g, 

. . . , z^^, . . . , z%j = uiz^ - 4)^-1 1[{4 - z^r' m< - (e-s) 

i<j k<l rs 

For functions on the sphere, the number of fiux seen by e and by o particles must be the 



same, but this is not the case for the function ( |6.3| ) as it stands, unless A'",, is zero. If A^^o > 



it is necessary to multiply the function by additional factors, and the space of these factors 
may be parametrised by viewing them as quasihole factors U„{w'') = Hii^f — w'^), where 
cr = e or o, acting on particles of either spin. On multiplying in Uf. factors of f/g, rio factors 
of Uoi one finds for the fiux seen by e and o respectively, 

N^ = {q-l){N,-l) + qN, + n, (6.4) 
= {q + l){No-l) + qN, + no (6.5) 

since the fiux must be the same for both; the second line, which is the fiux seen by the o 
particles, applies only if A'o > 0. 

We wish to analyse the situation A^^ > q{N — 1) — 1, that corresponds to the Pfaffian 
ground state or the same plus quasiholes. If we write 

^N^ = N^-[q{N (6.6) 

which was denoted n in Sec. |I|, then we find the g-independent equations 

ne = N,-2 + AN^, (6.7) 
no + No = AN^, (6.8) 

where again the second equation does not apply if A'o = 0. Now we see that if AA",^ < 0, 
we must have Aq = since rio > 0. So in this region, all zero-energy states contain only e 
particles. In particular, this includes the Pfaffian state and all states degenerate with it at 
AN^ = 0; for g = 1, this space of states includes all states in which all the bosons have spin 
e, and the same is true for the fermions at g = 2, as remarked above, since antisymmetry 
requires all states to include the factor Y[i<j{Zi ~ ^j) as a factor. The nondegenerate ground 
state of the model occurs at A",^ = (<? — 1)(^ ~ 1) (i-6- at AA",^ = — (A^ — 1) -|- 1), and for 
g = 2 is the LLL filled with e particles. For AA^^ > 0, the maximum number of o particles 
possible in the zero energy states is No < AN^. The numbers of zero-energy states, and their 
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angular-momentum decomposition, for each A^, N^p for this model can now be obtained, as 
in other cases (see especially the case of the Laughlin state, and the 331 state in the bosonic 
language). 

Numerical study of Ho's model confirms the above discussion. Figures ^ and |10] show 
representative spectra with no quasiholes, for various values of 9, including the 331 point 
^ = 0, and for = 6, g = 2. We chose Vq = 1, = 0.5. For d not at the Pfafiian 
point, the zero-energy ground state at = q{N — 1) — 1, is nondegenerate. As the Pfaffian 
point is approached, a set of states approaches zero energy, and at the Pfaffian point, the 
degeneracies of the zero-energy states are exactly those of the mm'n system that we have 
just analyzed. Already at 9 = vr/S, one can identify most of the states that reach zero 
energy at the Pfaffian point. Calculations not shown in the Figures also confirm that states 
with quasiholes have the degeneracies, and the angular momenta, of those for the 331 state. 
We notice that, in addition to the zero-energy states at the Pfaffian point, which are fully 
explained by the above analysis, there are also some very low energy excited states, for which 
we do not at present have a detailed explanation. 

There is a simple way to remove the "excess" degeneracy at the Pfaffian point, without 
destroying the 331-like behavior of zero-energy states elsewhere. The Pfaffian, with all 
particles e, would be selected by the three-body Hamiltonian used earlier in the spinless case, 
if it acted on the e particles. We also observe, from the structure of the wave functions (p. 21), 
that in these wave functions no three particles have total angular momentum 3A"^/2— 3(g— 1). 
Therefore, the three-body operator which projects each set of three particles onto angular 
momentum 3N^/2 — 3(g — 1) and spin 3/2 (or its analog on the torus) annihilates the Ho 
states, and the quasihole states and ground states on the torus derived from it for all d 
(the quantum numbers are most easily derived by considering spin-1/2 bosons for g = 1; in 
these Ho states, three bosons are never found at the same point — if they were they would 
be in a symmetric spin state). If it is added to the Ho Hamiltonian, all the degeneracies of 
zero-energy states will be maintained, away from the Pfaffian point. At the Pfaffian point, 
the total Hamiltonian now selects, at A^^ = g(A^ — 1) — 1, the Pfaffian ground state as the 
unique zero-energy state, in which all spins are e. For smaller A",^ there are no zero-energy 
states, and for larger A^^ the zero energy states are just those of the Ho form ( |6^ ) and its 
generalizations. In these states, the paired fermions are all e, but the unpaired ones can 
be either e or o. Consequently, the degeneracies are again those of the 331 state, not those 
of the Pfaffian. This was of course inevitable by continuity, given that no states now come 
down to zero energy at the Pfaffian point. Numerical spectra confirm these predictions, 
as shown in Fig. ^ for quasiholes, and in Fig. |T2| for 4 quasiholes. In these figures, the 
coefficient of the three-body projection operator is 1. 

We have now arrived at a Hamiltonian which has the degeneracies of the 331 state at 
all d vectors. Yet the Ho model was supposed to represent the effect of tunneling, which 
should raise the o particles to high energy, and it was expected that the ground state for 
strong tunneling would be the Pfaffian. While the ground state we find at the Pfaffian 
point is the Pfaffian, the problem is that unpaired particles of either spin can be in the zero 
modes. But this can now be easily (and exactly) remedied by adding to the Hamiltonian 
the tunneling term —tJ2iCx,i- Since the tunneling term is diagonal in the e-o basis we are 
using, it simply splits the states we have found. Clearly, all states containing only e particles 
now have energy —Nt, and these have the degeneracies of the Pfaffian state for any number 
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of quasiholes, or for edge states, or on the torus. States containing o particles are higher 
in energy. If we consider the same model, as a function of d, and let t vary with d, such 
that t is zero at the 331 point and of order one at the Pfaffian point, then we can say that 
its lowest-energy states are known exactly at the 331 and Pfaffian points, but unfortunately 
not in-between. There must be one or more phase transitions as the Hamiltonian is varied 
between these limits. Notice that there is a surpising effect at the Pfaffian point, which has 
edge states including unpaired fermions of either spin when t is zero, but the o states obtain 
a gap when t > 0. 

When studying more general Hamiltonians that include tunneling, Ho's state would in 
general be a better choice of trial state than the original 331 state, and d can be used as a 
variational parameter. We wish to emphasize that a large overlap of the ground state of a 
Hamiltonian with a particular trial state, say, the 331 state (as in []18|) does not prove that 
the state is in the corresponding (say, 331) universality class; to show that, its properties 
must be calculated, and the thermodynamic limit taken, and this is a more difficult task 
numerically. Conversely, a low overlap would not prove it is not in that class. 

Our discussion still leaves the question of whether for large enough t there is a transition 
to the Pfaffian universality class, or to some other class. Because of the similarity in ground 
state structure of the 331 and Pfaffian states emphasized in Ref. P, p!3| , pil| , a second order 
transition seems to be a possibility. Another possibility is that the Halperin class of paired 
states is involved, those for which a trial state can be constructed by first pairing the particles 
into charge-2 bosons, then forming auh = 1/Aq = 1/8 state of the bosons p|,^. This is clearly 
another abelian quantum Hall state, for which the edge excitations would consist only of 
the U(l) density-fluctuation chiral scalar boson modes, in different charge sectors; it would 
lack the gapless fermion excitations characteristic of the paired states discussed here [|]. 

While we cannot rule out the Halperin type of paired state, and all of these states 
might describe the universality classes of various Hamiltonians, even within the LLL, we 
can point out an experimental signature that will distinguish the former from the others. 
First we recall that, for the Laughlin states, the Luttinger liquid at the edge leads to the 



local density of states N{uj) ~ li^l^""*^ for fllling factor 1/q (see Ref. [31| for a review). The 
corresponding tunneling differential conductance for a point contact at zero temperature 
will be dl/dV ~ V^'^. The exponents are simply related to the scaling dimension A 
of the operator that creates a single electron in the low energy theory of the edge, as 
dl/dV ~ and also to the statistics 9 of the excitation by 9/27i = A (mod 1). In the 

theory for the Laughlin states, the charge n/q operators at a single edge have dimensions 
n^/2g (n must be an integer); the non-integral charge operators must be combined with 
similar operators at the other edge (see Refs. or |^). Thus A = g/2 for the electron. In 
the charge-2 boson system, the charged operators that create charges at the edge within the 
low-energy theory (subject to some straightforward conditions discussed in Ref. [§]) have 
scaling dimension ra^/Sg when they create charge ±n/2g (in electron units) at the edge at 
fllling factor 1/q, again for n integral. As in the old argument of Tao and Wu in either 
theory the operator of charge 1 has 6 = 27rA = vrg, which is Bose statistics for q even, and 
cannot represent an electron (similarly, for q odd, it cannot represent a charge-1 boson). In 
addition, such an operator cannot be used at a single edge in the charge-2 boson theory, 
but must be combined with another charged operator at the other edge, or with some other 
operator at the same edge. Only operators with even-integer charge can act on one edge. 
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Therefore, at low energies, (that is, bias voltages), tunneling into the edge from outside will 
be impossible at these filling factors v = 1/q for this charge-2 boson universality class. To 
make an electron on one edge, the charged operator must be combined with an operator 
making an unpaired fermion (or BCS quasiparticle) , which in the present case would exist as 
excitations, but would have a finite energy gap, even at the edge; thus the tunneling current 
will be zero below a threshold voltage. On the other hand, in the paired states like the 
Pfafiian, HR and 331 states, the fermion excitations are gapless at the edge, and a nonzero 
tunneling current with a power-law dependence at small bias voltage should be observable 
(the power law can be calculated from the theories in Ref . ; A has an extra contribution 
from the fermions, to yield V^, for both the Pfaffian and 331 cases, and V^^^ for HR). These 
arguments are for leading order in the tunneling. It is possible that some sort of higher-order 
tunneling process could transfer two electrons into the edge, through a virtual transition to 
one or more higher-energy states; as this does not require any fermions to be created except 
virtually, this would give a current at arbitrarily low bias voltage. However, the exponent 
would be related to the scaling dimension (and statistics) for the charge-2 operator in the 
low-energy edge field theory of the charge-2 boson state, and so the power law would be 
V'^'^~^\ thus for q = 2 the current would be much lower than for Pfafiian, 331 and HR states, 
and still clearly distinguishable. For another discussion of the experimental consequences 
of the Halperin-type charge-2 boson Laughlin state, in connection with even-integer filling 
factors, see Ref. |^; our formulas for the charge-2 boson universality class apply for the 



conductance in this case also, with q = 1/2, to give dl/dV ~ V. 



VII. CONCLUSION 



In summary, we have obtained a full description of the quasihole states of several paired 
FQHE states, for the Hamiltonians for which the exact ground states are known. The 
degeneracies found in the Pfaffian and HR cases are as required for nonabelian statistics. 
For the 331 states, the statistics are abelian, and the degeneracies are due to the layer 
index. Ho's model was found to be pathological at its Pfaffian point, but the pathology 
was removed by adding a three-body term to Ho's Hamiltonian. With tunneling added 
also, the Pfaffian state was recovered, but the model was no longer exactly soluble for 
the low-energy states at intermediate points in the parameter space between the 331 and 
Pfaffian points. In Appendix p, the permanent state was also considered, which is another 
candidate for nonabehan statistics, but should probably be rejected because of its proximity 
to a ferromagnetic ground state, and its correspondingly gapless nature. It remains to 
be seen whether the other nonabelian states, though not apparently close to an obviously 
gapless state, are in fact stable against small generic perturbations. This is an important 
outstanding issue to which we hope to return elsewhere. It probably requires an analytical, 
field-theoretic technique to settle it in general, which should be a theory that describes the 
paired condensate, and not just a Chern-Simons theory of the low-energy sector containing 
nonabelian statistics. In the meantime, we have pointed out in Sec. [V^ how the different 
paired states, the Halperin-type state of charge-2 bosons, and the Pfaffian and HR types with 
nonabelian statistics, can be distinguished in a point-contact tunneling experiment. As for 
an actual demonstration that adiabatic transport of quasiparticles does produce nonabelian 
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statistics in some systems, that also will have to be left for treatment elsewhere. 
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APPENDIX A: ELLIPTIC FUNCTIONS 



In general, theta functions with characteristics are defined as 
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where the n sum is over all the integers and a and b are real. From the definition we obtain: 



{z + 1\t) = e'^'^^'d 
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Consequently we can restrict a, b to lie between and 1. 

The standard Jacobi theta functions |^ are particular cases of those above. There are 



four of them: '(9i(2;|r) = 
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and {Id = 
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is of particular importance, since it is odd under z — z, and so has a zero at the origin. 

a 



. Of these. 



In general, all ^ 



can be related to each other by shifts of z. 



The identity used to rewrite the Fcm found in Ref. |2^ in the form given in Sec. |V| can 
be obtained, by shifts of z, from the simplest version 
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Y[ i^siz - r/M\T) = i^s{Mz + (M - l)/2|Mr)r/*^(r)/r/(Mr), 



(A3) 



r=l 



which is obtained by writing the theta functions on the left-hand side in the product form 
2| 



^^(z\t) = n (1 - e'""")(l + 2e2"'"(2"-') cos27r2 + e^'^^^^^n-i) 



(A4) 



n=l 



and doing the r product first. Here rj{T) is the Dedekind rj function. 



n=l 



(A5) 
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The Weierstrass elliptic function can be defined by p9 



Piz\T) = ^ + \- ^ — (A6) 



where the prime indicates that m = n = is to be omitted from the sum. It can be shown 
29| that p is invariant under + z^z + t. To obtain functions that have a 



double pole at the origin, like the Weierstrass function, but are antiperiodic, we may use 
the definitions 

PaAAr) = ^(-1)™^+"^- i (A7) 

^ [z-m- nry 

where a, b are integers. depends on a and b only modulo 2. For a, b both even, the sum 
is not convergent, which is why p was not defined this way. For the other three cases, the 
series converges, and it is clear that the functions have the periodicity properties that we 
denoted H — , — h, in the text. 



APPENDIX B: PERMANENT STATE 

In this Appendix, we will introduce a Hamiltonian for which a certain state containing a 
permanent (other such states can be found in Ref. 0), is the unique zero-energy eigenstate 
of maximum density, and summarise results for the quasihole and edge excitations, and for 
the ground states on the torus. We also describe a relation with fully spin-polarized states 
and their skyrmion excitations, and argue that the permanent state is at a phase transition 
from ferromagnet to paramagnet. 

The permanent state is a spin-singlet ground state of spin-1/2 fermions for q odd, and 
of spin-1/2 bosons for q even. It can be viewed as spin-singlet p-wave pairing of composite 
bosons [^]. The Hamiltonian for the simplest case, g = 1, is a three-body Hamiltonian, 
which penalises the closest possible approach of three spin-1/2 fermions. On the sphere, 
three particles are at their closest, consistent with Fermi statistics, when the total orbital 
angular momentum for the three is 3N^/2 — 1, and the total spin is 1/2. Our three-body 
Hamiltonian is therefore taken to be a positive number times the projection operator onto 
this multiplet of states for three particles, summed over all triples of particles: 

^ = E VP,,k{^N^/2 - 1, 1/2). (Bl) 

i<j<k 

We have verified numerically that this does produce a unique many-particle state at zero 
energy at the value that corresponds to the permanent state a.t u = 1. As the permanent 
state is a zero-energy state for this Hamiltonian by inspection (and is non-zero), it must be 
the state obtained numerically. As for the Pfaffian state in Sec. |T|, a suitable Hamiltonian, 
consisting of a combination of spin-independent two-body projection operators onto angular 
momenta N^, — 1, . . . , A^",^ — g + 2, and a certain three-body projection onto angular 
momentum 3A'<^/2— 3(g— 1) — 1 and spin 1/2, can be constructed for which the generalizations 
of these states to g > 1 are again the complete set of zero energy states. 
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In our much-abused notation, the permanent state is defined by the wave function 



The fact that it represents a singlet is most easily seen by viewing it as singlet pairs of 
composite bosons of spin 1/2. It is totally antisymmetric for q odd, symmetric for q even. It 
resembles the 331 state, but with the sign of the permutation omitted. Such a summation 
over permutations defines the permanent of a matrix, perm M = J2a Hi , in which M 

is an L X L matrix with elements Mij, and the sum is over all members a of Sl- 
States with 2n quasiholes can be written down in the now-familiar manner: 



B-t B\ 



E n(4(.))"^-n(^ m)™' n 



G5jv^ fc=l 1=1 r=l K^ai^B^+r) ^p(B^+r)) 

x\{{z,-z^)\ (B3) 

in which $ is as in Sees. II-IV, and n^, m/ must be < n — 1. In these states the fiux is 
A*"^ = q{N — 1) — 1 + 72. It is clear that these are all zero-energy states for the three-body 
interaction eq. (|B1| ) at g = 1, and its generalization to g > 1. The counting of states is most 
similar to the 331 state, but because the unpaired particles are here spin-1/2 bosons, the 
Pauli principle restriction on the number of unpaired particles does not apply, and there is 
no upper limit, except that the numer of unpaired particles cannot exceed the total number 
of particles. For n > 1, the number of states for fixed if 's, and a fixed number of broken 
pairs, is that for B = + Bj^ unpaired bosons, in 2n orbitals (including spin f or |), which 
yields 




1 



(B4) 



Including the positional degeneracy of the quasiholes, and summing over B as for the other 
states, gives 




{N - B)/2 + 2n' 
2n 



(B5) 



The sum over B satisfying (—1)^ = (—1)^ diverges as ^ oo. Note that for > 1, it 
is possible to break all the pairs and make all the spins polarized, still with zero energy. 
In other words, the quasihole states of the spin-polarized u = 1 state form a subset of the 
zero-energy states for the three-body Hamiltonian. 

As we have seen in the main text, the edge states are closely related to the bulk quasihole 
states. For the permanent state, we will just state that there are 4g sectors of edge states 
involving unpaired bosons, much like the other examples in Ref. 0. In the twisted sector, 
there are zero modes which can be occupied with arbitrarily many bosons of either spin. 
Turning to the ground states on the torus, there are again 4g sectors for A^ even. The 
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ground states for k 7^ are an obvious generalization of those for the 331 state, containing 
a permanent instead of a determinant. These states are again spin singlets. For k = 0, the 
construction that worked for 331 made essential use of antisymmetrization, and does not 
work here. Instead, the only possibility is to break all the pairs and put the bosons in the 
constant single-particle state on the torus. The spin state is then totally symmetric, so we 
obtain a spin N/2 multiplet of states. This construction also works for odd. 

The large degeneracies of states in certain sectors in this system make sense in the 
interpretation in terms of correlators in CFT. The theory relevant to the permanent is a 
/3-7 ghost system [Q, where (3 and 7 are free bosonic fields of conformal weight 1/2, so that 
this theory is nonunitary. It realises Kac- Moody symmetry at level = — 1/2 (not to be 
confused with the vector k on the torus) . The latter theory includes the "spin" fields for the 
/3-7 system, which behave as infinite-dimensional fractional-spin representations of SU(1,1) 
symmetry, that is related to, but not quite the same as, the SU(2) symmetry in which we 
are interested. Consequently the infinite degeneracies when quasiholes are present, or in the 
ground state on the torus in the corresponding sector, are not surprising. 

Another interesting question is the excitation spectrum when the zero-energy states of 
our Hamiltonian are fully (or even just macroscopically) spin-polarized, as occurs in the 
presence of 2 quasiholes or on the torus. Since the ground state breaks the spin-rotation 
symmetry, and the Hamiltonian is short-range, we expect low-energy spin-wave excitations 
to exist, in which one or more spins are flipped. For a generic Hamiltonian, one would expect 
these to occur at low but nonzero energy, with a gapless quadratic dispersion relation as in 
a quantum ferromagnet. Then the system would not be fully gapped for spin excitations, 
unlike (we believe) the other systems studied in this paper. It is interesting to note the 



relation with the u = 1 spin-polarized system that has been much studied recently [34 



For a two-body Hamiltonian consisting of a projection operator onto zero relative angular 
momentum for two particles of total spin 0, the spin-polarized filled Landau level state, 
and quasihole states including reversed spins (skyrmions) are zero-energy eigenstates (note 
that the number of quasiholes for this system is identified with the number of flux added 
to the polarised ground state, which is smaller by 1 than the number of flux added to the 
permanent ground state). For each number of added flux, these are exactly the same as the 
states above with no unbroken pairs. There are also the expected spin-wave excitations at 
nonzero energy. However, in the case of the quasihole states of the Hamiltonian we have 
studied here, we already know that there are other states with some additional reversed 
spins at zero energy; they are simply the states where not all the pairs are broken. Clearly, 
these states would not be zero-energy states for the two-body Hamiltonian. It is tempting 
to identify them with a subset of the spin-wave states. 

From the summand in (p5|), we can obtain the number of zero-energy states for each 
number of unbroken pairs, {N — B)/2. It is instructive to begin with the case n = 1 
that corresponds to the sector containing the fully-polarized filled Landau level state. As 
we increase (A^ — B)/2 from zero, we expect that the total spin must decrease due to the 
formation of singlet pairs. The first binomial coefficient in eq. ( [B5| ) is equal to B + 1 which is 
the degeneracy of a single multiplet of spin S = B/2. The second binomial coefficient in eq. 
( [B5| ) is the orbital degeneracy of the quasiholes of the permanent, which we here interpret 
as the number of ways of placing (A^ — B)/2 bosons in 2n + 1 = 3 orbitals. The bosons 
are to be viewed as the spin waves. The three orbitals form an L = 1 multiplet. The spin 
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wave excitations in general could have angular momenta 0, 1, . . . ; the L = ones simply 
represent a global rotation of the spin and have already been counted in the degeneracy of 
each spin multiplet. We conclude that, for our Hamiltonian, L = 1 spin waves have zero 
energy at n = 1. The other spin wave states would have to be obtained by a collective 
excitation of the condensate of singlet pairs, or spin waves (depending on our point of view), 
which excites one or more of them to higher angular momentum. A similar picture holds for 
n > 1. Based on these arguments, we do expect a gap less branch of low-energy spin waves 
to be present in the spectrum of elementary excitations of the spin-polarized zero-energy 
states for our Hamiltonian, on the sphere, and also on the torus. 

The three-body Hamiltonian ( [BID above has been diagonalized numerically for q = 1 and 
N up to 12, with V = 1/18, for various numbers of quasiholes; results for n = and n = 1 are 
shown in Figures |T3| and |T4|. For n = 1, the degeneracies of the zero-energy quasihole states 
have been confirmed, and in addition (see Fig. |l^) there are low-lying states, so the system 
is not obviously gapped. In fact, by examining the states at n = 1 with = N/2 — 1, that 
is one less than the maximum value, which we expect to be the single-spin-wave states of the 
ferro magnet, we obtain a dispersion relation of the spin waves, which is shown as the lower- 
right inset in Fig. |1^ It has the expected form of a finite-size version of a gapless branch of 
states, and has zero energy for both L = and 1; the latter property may mean that the 
dispersion relation has vanishing coefficient of wavevector-squared in the thermodynamic 
limit (this coefficient is proportional to the spin stiffness in the ferromagnet). We note that 
these states penetrate quite far into the full spectrum. The apparent slight gap above the 
zero-energy states for L 7^ 5 is in fact just a finite size effect, since at least the Sz = N/2 — 1 
states must form a gapless branch as ^ 00. The low- lying states in the full spectrum 
should be other multi-spin-wave states. 

For the ground state sector of the permanent, there is a nondegenerate zero-energy 
state, as in the other cases studied in this paper, but there is also an apparent gapless 
branch of states at low L (see Fig. O). From the point of view of the ferromagnet, this N^ 
value represents a quasielectron state, which would be an antiskyrmion. For the two-body 
Hamiltonian, the antiskyrmions form a set of states with L = S, but which do not have 
zero energy or exactly soluble wave functions. The lowest-energy states in Fig. |T3| have 
L = S for L = 0, 1, 2, 3, but S = 1 for L = 4, though here an 5* = 4 state lies at slightly 
higher energy. Moreover, the L = ground state has overlap-squared 0.81 with that for the 
Coulomb interaction at the same A^, A"^. There is thus some evidence that this branch of 
states represents something similar to an antiskyrmion. 

The proximity to the ferromagnet, albeit at a different number of flux (A^^ = g(A^ — 1) for 
the ferromagnet), suggests that the system is at the transition to the ferromagnet. Indeed, 
if the two-body Hamiltonian containing Vq only (for g = 1, or the usual generalizations 
for other q) is added to the three-body, then for Vq > the ferromagnet and its quasihole 
excitations (skyrmions) will be the only remaining zero-energy states, and the L = 1 spin 
waves (at A"^ = g(A^ — 1)) will have non-zero energy. On the torus, the 4g-fold degenerate 
sectors will be split to leave only the q states with k = 0, which have maximum spin A^/2. 
For Vq negative, we expect the splittings to reverse sign, and the ground states on the torus, 
or on the sphere with quasiholes will presumably be unpolarized; in this region no exact 
wave functions are available for the lowest energy eigenstates. In other words, the polarized 
and unpolarized states will differ in energy density. We expect that the unpolarized states 
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are paired, and indeed the attractive pseudopotential Vo < should favour the pairs. In 
view of the higher energy density of the polarized ground state on the torus, we guess that 
these states arc no longer in the same universality class as the permanent state, but may be 
a simple Halperin-type state, a Laughlin state of charge-2 boson pairs. This probably would 
occur because the attractive potential decreases the size of the pairs, and modifies the pairing 
function from the simple form (zj — Zj)~^ found hitherto. As nonabelian statistics probably 
relies upon the long-range character of this part of the wave function, it could disappear 
under this perturbation. In any case, the ferromagnetic order parameter is constant for 
I/q > 0, and will vanish for Vq < 0, which indicates a first-order phase transition. 

The transition also has a simple interpretation in terms of composite bosons. In the 
permanent state, the bosons are paired, but when one flux quantum is added, there are 
broken pair states of the same energy. In these states, the unpaired composite bosons 
occupy the L = zero mode, and can be viewed as a Bose condensate. Since they carry 
spin 1/2, such a condensate is necessarily a ferromagnet, and when more flux are added, the 
skyrmion zero-energy states are obtained. Thus the permanent three-body Hamiltonian can 
have pairs but the bosons can also unpair and form a condensate. In Bose liquid systems, the 
condensation of single bosons is the usual occurence. It seems that the permanent state is on 
the borderline between a single-particle Bose condensate and a condensate of pairs only. It 
is possible that, while the ferromagnetic order indicates a flrst order transition, the behavior 
of the pair order parameter (speciflcally, the size of the pairs) , or of some properties on the 
ferromagnetic side, possibly related to a spin stiffness going to zero, could be characteristic 
of a second-order transition, with a length that diverges at the transition point. 
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FIGURES 

FIG. 1. Two-particle correlation function g{r) for the Pfaffian state, with = 2{N — 1) — 1 
(i.e. V = 1/2), for = 10 (dashed line), 12 (dot-dashed line), 14 (solid line), versus great circle 
distance on the sphere. The curves for = 12 and 14 are almost indistinguishable. 

FIG. 2. Spectrum of the three-body Hamiltonian for the Pfaffian state of fermions with v = \ /2, 
for = 12, Nm = 21, that is, no quasiholes. The inset enlarges the low-lying levels. 



FIG. 3. As in Fig. |, but with N = 10, = 19, that is, 4 quasiholes. 

FIG. 4. Low-lying excited states of the three-body Hamiltonian for the Pfaffian ground state 
at = 1/2 (iV0 = 2(A^ - 1) - 1) for iV = 10 and 12, plotted against k = L/R. 

FIG. 5. As in Fig. g, but for the energy A of the lowest-energy excited state versus 1/A^. 

FIG. 6. Ground state energy of the three-body Hamiltonian for the Pfaffian state at 
A''^ = 2(A^ — 1) — 2, that is, v = 1/2 with two quasielectrons added, plotted against 1/A^. 

FIG. 7. Spectrum of the hollow-core model for the HR state of fermions with i/ = 1/2, for 
N = 8, = 12, that is, no quasiholes. The inset enlarges the low-lying levels. 

FIG. 8. As in Fig. 0, but with A^ = 6, A''^ = 12, that is, 8 quasiholes. 

FIG. 9. Spectrum of the Ho model, for four values of 6 which parametrizes the d vector, all for 
A^ = 6, A'^^ = 9, that \s v = 1/2 and no quasiholes. 

FIG. 10. Same as Fig. |^, but enlarged to show low-lying levels. 

FIG. 11. Spectrum of the Ho model plus the three-body interaction at the Pfaffian point 
6 = vr/4, again for A^ = 6, A'^^ = 9. The inset enlarges the low- lying levels. 



FIG. 12. As in Fig. 11, but A^ = 6, A^<i = 11, that is, 4 quasiholes. 



FIG. 13. Spectrum of the three-body interaction for the permanent state of fermions with 
V = 1, for A^ = 12, A'^^ = 10, that is, no quasiholes. The inset enlarges the low-lying levels. 
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FIG. 14. As in Fig. |T^, but with N = 10, = 9, that is, 2 quasiholes. Upper-left inset: 
low-lying levels. Lower-right inset shows only the states with Sz = N/2 — 1, i.e. single spin flips 
(or spin waves) of the fully-polarized state. 
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TABLES 



TABLE I. Numbers of multiplets of states of total angular momentum L at zero energy for the 
three-body Hamiltonian for the q = 2 PfafHan state on the sphere, for = 2(A — 1) — 1 + n, i.e. 
2n quasiholes. 

TABLE II. Numbers of multiplets of states of total angular momentum L and total spin 
S at zero energy for the hollow-core Hamiltonian for the q = 2 HR state on the sphere, for 
= 2{N — 1) — 2 + n, i.e. 2n quasiholes. 
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